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POSSIBLE NEW EXPERIMENTAL TEST OF GENERAL RELATIVITY THEORY* 


L. I. Schiff 
Institute of Theoretical Physics, Department of Physics, Stanford University, Stanford, California 
(Received February 11, 1960) 


In a paper now in process of publication,’ it is 
argued that only the planetary orbit precession 
provides real support for the full structure of 
the general theory of relativity. The other two 
of the three “crucial tests,” the gravitational red 
shift and deflection of light, can be inferred cor- 
rectly from the equivalence principle and the 
special theory of relativity, both of which are 
well established by other experimental evidence. 
It is also pointed out that a terrestrial or satellite 
experiment that would really test general rela- 
tivity theory would have either to use particles 
of finite rest mass in such a way that the equa- 
tion of motion can be confirmed beyond the New- 
tonian approximation, or to verify the second- 
order deviations of the metric tensor from its 
Minkowski form. 

In an attempt to devise a feasible experiment 
that might accomplish one of these objectives, 
we have calculated the properties of a spinning 
test particle (torque-free gyroscope). We start 
from the covariant equations of Papapetrou’ for 
the motion of the center of mass and the spin 
angular momentum, generalized by inclusion of 
anongravitational constraining force F, and work 
to lowest order. The motion of the center of 
mass in the gravitational field of the rotating 
tarth is then described by the Newtonian equa- 
tion 


m(d¥ /dt) = -(GmM/r)r +F, (1) 


where m is the rest mass of the particle, fr is 
its coordinate, V =dr /dt is its velocity, G is the 
Newtonian gravitational constant, and M is the 
mass of the earth. The spin angular momentum 








vector measured by a co-moving observer, 8°, 
obeys the equation 


dS° /dt = 2 x8°, (2) 
where 
& = (F xv) /2me? + (3GM/2c?r)(F xv) 
+(GI/c?r9)[3(G-F)F/r?-5); (3) 


I=2MR?/5 is the moment of inertia of the earth 
of radius R, assumed to be homogeneous, and 

@ is its angular velocity vector. The first term 
on the right side of (3) is the Thomas preces- 
sion,® which is a special relativity effect. The 
other two are the lowest order effects of general 
relativity; the second term arises whether or 
not the earth is rotating, and the third term is 
the earth rotation effect of Lense and Thirring.* 
While the second term involves the first-order 
deviations of the metric tensor from its Minkow- 
ski form, which can be calculated without the 
use of general relativity,’ it also depends on the 
equation of motion of matter of finite rest mass 
beyond the Newtonian approximation. It is there- 
fore a genuine consequence of general relativity. 
The same is true of the third term, since in 
addition it depends on off-diagonal space-time 
components of the metric tensor. 

Equations (2) and (3) may be obtained either 
from the standard or the isotropic form of the 
Schwarzschild line element, and using for the 
supplementary condition on the angular momen- 
tum tensor either that of Corinaldesi and Papa- 
petrou’ or of Pirani.* The equation of motion of 
the spin in the nonrotating, earth-centered co- 
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ordinate system looks quite different in these 
four cases, but they all agree when expressed 
in terms of the spin measured by a co-moving 
observer. It should also be remarked that the 
corrections to Eq. (1) that arise from the spin 
are unobservably small in any realizable situa- 
tion. 

It follows at once from the form of Eq. (2) that 
the magnitude of the spin angular momentum 
measured by a co-moving observer is constant 
in time. Thus if the moment of inertia of the 
spinning particle does not change, the angular 
velocity of rotation is constant, and the spinning 
particle behaves like a clock which can be set to 
any desired frequency. This frequency exhibits 
Doppler and gravitational shifts when observed 
from outside, just like that of a more conventional 
clock. It is possible that its frequency stability 
could be made to compare favorably with those 
of other types of precision clocks. It also fol- 
lows from (2) that a number of spinning particles 
with various magnitudes and directions for their 
angular momentum vectors maintain fixed angles 
of these vectors with respect to each other. The 
vector 2, which in general is not constant, is 
their common angular velocity of precession with 
respect to the external “fixed stars”; in com- 
paring their directions with the outside world, 

a correction must of course be made for aber- 
ration whenever V #0. 

If a spinning particle is in free fall, as ina 
satellite, then F=0. For an orbit in the earth’s 
equatorial plane, for example, 


& = (3GM /2c*r)&, - (2MGR?/5c?r)5, (4) 


where @,=(f x¥)/r* is the instantaneous orbital 
angular velocity vector of the particle. The 
minus sign in Eq. (4) deserves some comment. 
The third term of Eq. (3) tends to cause a spin- 
ning particle to precess in the same direction as 
the rotating earth at the poles (F parallel or anti- 
parallel to @), but in the opposite direction at 
the equator (fF perpendicular to @). This is 
physically reasonable if we think of the moving 
earth as “dragging” the metric with it to some 
extent. At the poles, this tends to drag the spin 
around in the same direction as the rotation of 
the earth. But at the equator, since the gravita- 
tional field falls off with increasing 7, the side 
of the spinning particle nearest the earth is 
dragged more than the side away from the earth, 
so that the spin precesses in the opposite direc- 
tion. 

If the center of mass of the spinning particle is 
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constrained to remain at rest with respect to the 
rotating earth, as in an earth-bound laboratory,’ 
then 


¥V=@xr, dv/dt=@xv. (5) 


The required constraining force F can then be 
found from (1) and (5), and substituted into (3), 
When the particle is at the surface of the earth 
at latitude A, the precession angular velocity 
may be written in the form 


@ =[(4gR /5c?)(1 + cos*d) - (wR? /2c?)cos?a Ja 
+(4gsinA/5we?)(@x¥), (6) 


where g=GM/R? is the acceleration of gravity 
at the surface of the earth. Only the square 
bracket term in Eq. (6) gives rise to a secular 
precession of the spin axis, and the second part 
of it is very small compared to the first. Thus 
to good approximation, a particle with spin axis 
perpendicular to the earth’s axis precesses at 
the rate 27(4gR /5c”)(1 +cos?A) =3.5 x 107*(1 +cos*) 
radians per day. It also follows from Eq. (3) tha 
the corresponding effects caused by the sun and 
moon are negligivly small in comparison. 

A secular precession of 6 x10~° radian per day 
would be very difficult, but perhaps not impos- 
sible, to observe. Professor W. M. Fairbank 
and Professor W. A. Little of this department 
are exploring the possibility of using for this 
purpose a gyroscope that consists of a super- 
conducting sphere supported by a static magnetic 
field.* Such a gyroscope would also be of interest 
as a device for performing experiments in low- 
temperature physics. If it could be made to 
operate exceedingly well, it might in addition be 
used for an experimental test of Mach’s princi- 
ple, by comparing the orientation of its axis 
with a field of “fixed stars” over a period of a 
year or so. Most of the experimental difficulties 
that seem to arise with a high-precision gyro- 
scope are greatly reduced if the gyroscope does 
not have to be supported against gravity. This, 
together with the fact that w, is generally much 
larger than w, suggests that experiments of this 
type might be more easily performed in a satel- 
lite than in an earth-bound laboratory. 

A full account of this work will be submitted 
shortly for publication in the Proceedings of the 
National Academy of Sciences. 





“Supported in part by the U. S. Air Force through 
the Air Force Office of Scientific Research. 
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BUMPY TORUS* 


Gordon Gibson,! Willard C. Jordan,t and Eugene J. Lauer 
Lawrence Radiation Laboratory, University of California, Livermore, California 
(Received February 1, 1960) 


Magnetic field configurations resulting from 
various coil arrangements are being used for 
thermonuclear studies. In the absence of par- 
ticle interactions, charged particles which oc- 
cupy certain regions of velocity space are con- 
tained indefinitely provided the particles behave 
adiabatically. However, there are regions in 
velocity space in which particles are not con- 
tained; e.g., in the mirror machine if the velo- 
city vector of a charged particle makes a suffi- 
ciently small angle with the magnetic field the 
particle escapes through a mirror.’ Hence, the 
mean containment time is limited to approxi- 
mately the time, 7s, for small-angle scattering 
tocumulate to a large angle insofar as multiple 
scattering is the dominant loss mechanism. 

A uniformly wound torus would eliminate the 
end losses; however, another difficulty is intro- 
duced. Because of the drift of the charged par- 
ticles resulting from the gradient of the magnetic 
field, the particles rapidly impinge on the walls 
of the device. 

A simple case to illustrate the drift motion 
both ina magnetic mirror machine and in the 
torus is one where the velocity, v, of the par- 
ticle is normal to the magnetic field, H, and the 
longitudinal gradient of the field is zero. The 
magnetic moment, yu, of the particle motion is 
an adiabatic invariant?; therefore .=W,/H 

=;mv,*/H =constant, and if the particle’s en- 
ergy, W, is a constant the particle’s guiding 


center moves along a path defined by H = constant. 


The particle drifts are illustrated in Fig. 1. 

A rotational transform’ tends to reduce the 
loss due to drifts in the torus. However, the 
rotational transform may be ineffective if the 
particle motion along the field lines is slow (re- 
lative to the drift velocity), with the result that 
the containment time would again be limited to 








FIG. 1. 
cross section, median plane. 
torus cross section. 


Particle drift paths. 


(a) Mirror machine 
(b) Uniformly wound 


approximately the scattering time 7,. The elec- 
tric field and other phenomena associated with 

a dense plasma would also affect the containment 
time. 

In this paper the single-particle motion in the 
magnetic field created by a circular array of 
circular current loops is investigated (see 
Fig. 2). The extent to which this geometry would 
have the advantage of closed precessional sur- 
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FIG. 2. The bumpy torus geometry and field mag- 
nitudes. The symmetry plane, S, is the plane of the 
sheet of paper. 
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faces which would not intersect the walls [see 
Fig. 1(@a)], and would also eliminate end effects,* 
is investigated. Because of the complexity of a 
3-dimensional calculation, the results in this 
paper are limited to an investigation of the in- 
tersection of the precessional surfaces with the 
symmetry plane containing the torus. The 
single-particle model is of special interest 
during the plasma buildup stage in the high-en- 
ergy injection schemes.*»*® The stability of a 
plasma has not been investigated theoretically 
for this geometry; of course if a plasma would 
build up, the cooperative effects could be ob- 
served experimentally. 

It is assumed that the following quantities are 
constants of the motion: (1) the kinetic energy, 
W; (2) the magnetic moment, » =W,/H; and (3) 
the action integral, J={p,di, evaluated over a 
period (as defined later) of the longitudinal mo- 
tion, where p N is the component of the momen- 
tum parallel to the field and di is the element of 
path length parallel to the field.” 

The IBM 650 is being used to calculate (1) the 
magnitude of the magnetic field, (2) the direction 
of the magnetic field, and (3) the action integral 
for particular values of the magnetic moment 
and starting coordinates. 

The code is set up to calculate these quantities 
in the symmetry plane S (see Fig. 2) for N 
equally spaced circular current loops (current J) 
with a ratio of major radius divided by current 
loop radius, R/r. Because of the symmetry it 
is only necessary to investigate the region be- 
tween a median plane and the plane of an adjacent 
current loop. 

Taking account of the constancy of the kinetic 
energy and of the magnetic moment, the action 
integral may be written as 


J =mv 3-2 sin?) dl, 
fir 


where mv is the momentum, Jd, is the initial 
field magnitude, and 6 is the initial angle that the 
velocity makes with the field direction. For a 
given value of »., this integral is calculated as a 
function of the radial position of the starting 
point in the median plane, M. The integrals are 
stopped either at the turning point where v i has 
decreased to zero or at the plane of the adjacent 
coil if there is no turning point.® 

Figure 3 shows the results of the action inte- 
gral calculations for seven values of py and also 
some of the field magnitude calculations that 
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FIG. 3. Adiabatic invariants and field magnitudes 
for a bumpy torus with N=16 and R/r=5.44. 


have been obtained for the case of N=16 and 
R/r=5.44. All of the action integral curves 
have a minimum (for J #0). Choosing a particu- 
lar curve of 1 =constant, then for a chosen value 
of J the two values of p (inner and outer), where 
the precessional surface intercepts the plane of 
symmetry S, are uniquely determined. It is 
most unlikely that these surfaces are not closed 
outside of the S plane; however, a 3-dimensional 
calculation is required for a rigorous proof. 
For the particular case of precession in the 
median plane M, J=0 and the intercepts of the 
precessional surface are not defined. However, 
as shown above, the precession for this special 
case is constrained to a path along which H =con- 
stant. The plot of H vs p/r in the plane S has a 
maximum (this must be a maximum for displace- 
ments out of the plane S also). Thus, inner and 
outer intercepts for the precessional curve 
(which must be closed) with the plane of sym- 
metry S may be found from the curve of H vs p 
in the median plane. In the limit as yp approaches 
the value characteristic of a particle trapped in 
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the median plane (holding the kinetic energy con- 
stant), the intercepts predicted by the J curves 
appear to approach those predicted by the 
median-plane H curve as they should. 

In case UV) is nearly zero at a mirror (transi- 
tion from dashed to solid lines in Fig. 3), J does 
not necessarily remain constant in time (see 
footnote*). Such a particle has v, =v =constant 
and so follows a closed precessional curve for 
which H = constant in the plane of the coil (i.e., 
at @=0). 

Another special case is for , =0. In this in- 
stance J is proportional to the length of the field 
line (over the path of integration), and since J is 
shown to have a minimum value in the S plane 
this is certainly a minimum with respect to dis- 
placements out of the S plane also. Hence, these 
precessional surfaces are closed. 

Scattering can change the value of » suddenly 
so that a vertical transition occurs in Fig. 3 
from one of the ». =constant curves to another at 
anearly constant p. After the transition the 
particle is constrained to the new precessional 
surface associated with the new values of yp and 
J until another scattering occurs. Since the J 
vs p curves do not all have the same shape and 
since their minima do not all occur at a single 
value of p, scattering can result in a new type of 
diffusion, that is, the particles can diffuse onto 
larger and larger diameter precessional sur- 
faces. (This diffusion is independent of the orbit 
diameter.) It is planned to investigate different 
values of N and R/r. Increasing N and R/r pro- 
portionally should enable one to approach the 
degenerate case of a linear series of mirror 
machines and hence reduce the effect of this new 
type of diffusion. Initial results with N=4 x16 
or 64 and R/r =4 5.44 or 21.76 indicate that the 
minimum number of 90° changes in pitch angle 
of the velocity that can change the diameter of 
the precessional surfaces by an amount equal to 
the coil radius is approximately three for par- 
ticle positions such that the distance from the 
axis is>r. Hence it is expected that the mean 


time for this change is about 3*7,, or 97g. 
We are grateful for the use of the results of 

T. G. Northrop and E. Teller prior to publica- 

tion. This work was carried out with the support 

and encouragement of C. M. Van Atta. The nu- 

merical work was done by J. L. Brady, T. H. 

Haratani, and B. Williams. 





"This work was done under the auspices of the U. S. 
Atomic Energy Commission. 

T Visiting from Westinghouse Atomic Power Depart- 
ment, Pittsburgh, Pennsylvania. 

t Visiting from Bendix Aviation Corporation, Re- 
search Laboratory Division, Detroit, Michigan. 

'E. Fermi, Phys. Rev. 75, 1169 (1949). 

2H. Alfvén, Cosmical Electrodynamics (Oxford Uni- 
versity Press, New York, 1948), p. 21. 

3L. Spitzer Jr., Phys. Fluids 1, 253 (1958). 

‘B. B. Kadomtsev, Plasma Physics and the Problem 
of Controlled Thermonuclear Reactions (Akad. Nauk. 
USSR, 1958), Vol. III, p. 285. A recent translation of 
this Russian paper indicates that G. I. Budker inde- 
pendently has suggested the possibility of stabilizing 
drifts by means of a “corrugated field.” In this paper 
Kadomtsev analyzes field configurations resulting from 
linear series of mirror machines connected by uniform 
field “sleeves.” 

5A. Simon, Phys. Fluids 1, 495 (1958). 

6G. Gibson, W. A. S. Lamb, and E. J. Lauer, Phys. 
Rev. 114, 937 (1959). 

'T. G. Northrop and E. Teller have recently shown 
that the particle drifts are such as to conserve the ac- 
tion integral for a general mirror-type field: Univer- 
sity of California Radiation Laboratory Report UCRL- 
5615, 1959 (unpublished). This has also been proven 
by Kadomtsev (see reference 5). 

®The justification for this assumption is that it makes 
J approach the same value in two limiting cases (hold- 
ing the kinetic energy and the magnetic moment con- 
stant): (1) for a particle that is trapped between two 
adjacent mirrors as the turning point approaches the 
mirror, and (2) for an untrapped particle as vy at the 
mirror approaches zero. A more detailed analysis of 
this case should be made because the longitudinal 
period becomes very long as the above limit is ap- 
proached, and the action integral theorem is not proven 
for this case. 
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EVIDENCE FOR THE QUANTIZATION OF CIRCULATION IN LIQUID HELIUM n* 
C. E. Chase 


Lincoln Laboratory, Massachusetts Institute of Technology, Lexington, Massachusetts 
(Received January 29, 1960) 


According to the theory of Onsager’ and Feyn- 
man,” rotating liquid helium II should contain an 
array of vortex lines arranged parallel to the 
axis of rotation and approximately evenly distri- 
buted throughout the liquid; the circulation 
around these lines should be quantized in multi- 
ples of h/m, where m is the mass of the helium 
atom. The lowest energy state of this system 
would arise when each line had only one quantum 
of circulation associated with it; the number of 
lines per square centimeter crossing a plane 
normal to the axis of rotation would then be 
2wm/h. Although there is some experimental 
evidence in favor of this model,*»* the actual 
quantization of the circulation has not yet been 
established in an entirely convincing manner.® 
This Letter reports preliminary measurements 
of the effect of rotation upon the critical heat 
current W, in a channel of moderate dimensions. 
The most straightforward interpretation of the 
results of these experiments suggests that quan- 
tized circulation is indeed a property of the ro- 
tating state, although the possibility of some 
other explanation cannot be positively excluded 
at this time. 

The apparatus (Fig. 1) consists of a thin-walled 
stainless steel tube C, 0.262 cm in diameter by 
5.16 cm long, one end of which is attached to a 
chamber A containing a heater H and a resistance 
thermometer 7,; the other end opens into the 
helium bath, where a second thermometer 7, is 
located. The entire assembly is surrounded by 
a vacuum jacket D. Heat currents are generated 
by the heater H, and the resulting temperature 
gradient is measured by observing the difference 
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FIG. 1. The apparatus. 


between 7, and T, on a 33-cycle differential 
thermometer bridge. C is packed with 70 stain- 
less steel wires, 0.025 cm in diameter, pro- 
ducing a large number of small parallel channels 
of irregular cross section (these wires have been 
omitted from the figure for clarity). The appa- 
ratus can be rotated about the axis shown in the 
figure at any speed up to about 6 rad/sec by 
means of a variable-speed motor mounted out- 
side the cryostat. Critical heat currents were 
measured by an ingenious method proposed by 
Vinen.* This method depends upon the fact that, 
when a heat current is switched on in undisturbed 
liquid helium, mutual friction only appears after 
a time 7 of the order of several seconds, but 
that this time is greatly reduced if a heat cur- 
rent W, greater than the critical value is initially 
present. By measuring 7 as a function of W,, it 
is thus possible to determine W, with a precision 
limited only by the size of the increments in W, 
used in the experiment. In the present case the 
resulting uncertainty in W, was usually + 0.001 
watt/cm’, although it was occasionally twice as 
large. Further experimental details will be 
given in a later publication. 

In Fig. 2 the critical heat current at 1.41°K is 
plotted as a function of w”?. The various sym- 
bols refer to different runs; vertical lines 
through the points indicate the experimental un- 
certainty. Open symbols represent points taken 
upon increasing the angular velocity to the indi- 
cated values, whereas points represented by 
solid symbols were taken after w had been re- 
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FIG. 2. Critical heat current at 1.41°K as a function 


of the square root of the angular velocity. Open sym- 
bols—increasing w; solid symbols—decreasing «. 
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duced from a higher value. It should be noted 
that, as w is increased, W, drops in a series of 
fairly sharp steps equally spaced along the ab- 
scissa. There appears to be some hysteresis at 
the first step, as indicated by the arrows; no 
such phenomenon has yet been observed at higher 
angular velocities, but detailed measurements 
have not been carried out there. 

Vinen’s®»” studies of mutual friction led him to 
propose that a heat current generates turbulence 
in the superfluid component, and that this tur- 
pulence takes the form of a tangled mass of vor- 
tex line. He further suggested that the critical 
heat current results from a competition between 
the mechanisms responsible for generation of 
vortex line and its annihilation at the walls of the 
channel. If these ideas are correct, it is plau- 
sible that the presence of vortex lines produced 
by rotation would alter the balance between these 
mechanisms in favor of the growth of turbulence, 
and thus decrease W,. A reduction in W, would 
then correspond to an increase in the amount of 
vortex line present in the channel, and the fact 
that this reduction takes place abruptly at cer- 
tain angular velocities would indicate that the 
amount of vortex line present increases in dis- 
crete steps. This suggests that within a given 
range of angular velocities only certain configura- 
tions of vortex lines are stable, and that for rates 
of rotation smaller than that corresponding to the 
first step in Fig. 2 no vortex lines are present at 
all. It is unlikely that such a situation would 
arise unless the circulation were strictly quan- 
tized, for otherwise the circulation around each 
line could vary continuously so that the necessary 
boundary conditions were satisfied at any angu- 
lar velocity. It appears that the number of stable 
configurations is in fact limited by the require- 
ment that they must rotate with the same angular 
velocity as the apparatus, each line moving in 
the velocity field of its neighbors and of the 
images in the channel walls; otherwise the lines 
would be rapidly annihilated at the walls or 
driven out the ends of the channel. 

If the circulation around each line is h/m, the 
equilibrium separation of vortex lines in a large 
volume of liquid rotating at an angular velocity 


w is 1=(h/2wm)”?. At the angular velocity of the 
first step, 0.4+0.15 rad/sec, this separation is 
0.034 +0.005 cm, which is comparable with, 
although somewhat larger than, the average 
width of the interstices between the wires of 
which the channel is composed. The values of / 
at the second and third steps are 0.0174 0.001 
cm and 0.0114 0.0002 cm, respectively. Inas- 
much as these lengths are in the ratio 1:4:4, it 
is tempting to suppose that the steps correspond 
to the appearance of one, two, and three paral- 
lel rows of vortex lines along the channel. Un- 
fortunately, it does not appear that such simple 
configurations would be stable against small 
disturbances or would rotate with the required 
angular velocity, unless they are constrained by 
friction with the top and bottom walls of the 
channel. 

The interpretation of these measurements 
would be much clearer if a single channel of 
simple geometry were employed. Similar ex- 
periments are therefore being carried out in a 
rectangular channel of dimensions 0.041 x0.109 
cm. Preliminary results with this apparatus are 
in full agreement with those described above. 

The author is grateful to Dr. E. Maxwell, Dr. 
D. Douglass, and G. Myers for many helpful 
discussions. 
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PLASMA TRAPPING IN CUSPED GEOMETRIES* 


Harold Grad 
Institute of Mathematical Sciences, New York University, New York, New York 
(Received November 10, 1959) 


The question of the trapping of particles in 
cusped magnetic containment configurations has 
been recently examined.’ The purpose of this 
Letter is to make accessible and extend certain 
early results which bear on this problem.” The 
trapping problem is, clearly, very closely re- 
lated to the question of containment. The latter 
problem has been treated in some detail for the 
cusped geometry’ (high 6 =2u.,)/B*) and toa 
lesser extent for the picket fence* (low 8). The 
possibility of creating a cusped plasma by in- 
jection using a plasma gun has been described 
several times® and has been tried experimentally 
on a small scale with inconclusive results.® 

The mechanisms are different for trapping of 
a single particle (i.e., a low-8 stream) or a high- 
B cohesive plasma burst; we consider first the 
simpler single-particle mechanism (this has 
been termed impossible; see reference 1). 

We consider a two-dimensional cusped con- 
figuration as shown in Fig. 1. In this simplified 
model the magnetic field is assumed to be totally 
excluded and a particle orbit consists of straight- 
line segments joined by helical arcs. Sufficiently 
close to a cusp, the orbit has the adiabatic in- 
variant’ (see Fig. 1 for notation), 


u* =2mvd +1m?(v? +w")/eB. (1) 


The value of u* will be essentially constant 
while a particle is close to a cusp, but it will 
change its value when the particle crosses the 
device and approaches a second cusp. The sim- 
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FIG. 1. Two-dimensional cusped geometry. 
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plest containment theory* is based on the assump- 
tion of a complete loss of memory on traversing 
the apparatus. In other words, the distribution 
function of particles approaching a cusp is Max- 
wellian without any gap in a loss-cone. A re- 
finement of this theory takes into account a 
certain amount of persistence in two successive 
values of u*.? Turning to the injection problem, 
we see that a complete loss in memory of the 
value of u* yields complete trapping of all in- 
jected particles. They will enjoy a lifetime 
which is, on the average, the same as that of 
any particle in the plasma, whether recently 
injected or not. It is only those particles witha 
good memory which cross to the opposite cusp 
and escape after but a single transit. 

A simple computation of this persistence effect 
(based on a calculation of the last encounter made 
with the original cusp before crossing to the 
opposite side) shows that all particles are trapped 
which are injected within a “trapped zone” (see 
Fig. 1) of width 


V2 
20=2(1 -qim) (2) 


where \ = m(v? +w?)”* /eB is the Larmor radius 
and M =u/(v? +w?)” is the injection “Mach num- 
ber.” To be trapped, the guiding center of an 
approaching particle should be within a distance 
6 of the cusp axis. With A/L fixed, there are no 
trapped particles when M is sufficiently large. 
The maximum width of the trapped zone is twice 
the Larmor radius. For most efficient trapping, 
the beam width should not exceed the Larmor 
radius. This formula is an approximation which 
is valid when either M or L/,d is large. 

A similar (but much more intricate) analysis 
can be performed for the case of axial injection 
through a point cusp in a three-dimensional 
geometry (Fig. 2). The maximum size of the 
trapped zone is in this case a disk of area 70° 
~L; this area has the same order of magnitude 
as the total trapping area offered by the line- 
cusp which can be analyzed two-dimensionally as 
above. An appreciable fraction of those particles 
which approach the point cusp within the distance 
6 will be trapped provided that M =u/(v? +w?)” 
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FIG. 2. Three-dimensional cusped geometry. 


is smaller than some numerical factor times 
(L/a)”®. 

The problem of computing the trapping of a 
high-8 plasma burst which can alter the mag- 
netic field is much more difficult. A simple 
fluid-like theory can be developed for the smooth 
motion of a plasma in a magnetic channel.® If 
we are guided by the similar problem in ordinary 
fluid dynamics, we can expect complicated phe - 
nomena such as shocks (as in a diffuser with 
supersonic entrance velocity), jet detachment 
(since the magnetic walls curve away sharply), 
vortices in the dead space surrounding the jet. 
Peculiar to the magnetic channel problem, it 
has been shown that there is no stable transition 
between subsonic and supersonic regions.* Thus 
the entry of the plasma into the low-field region 
isa complex matter. At the far end, one cannot 
conclude (as is done in reference 1), merely 
from the fact that the magnetic pressure is 
higher than the energy density in the beam, that 
the beam will be turned back. For one thing, 
in the direction parallel to the magnetic field, 
the Maxwell stress is a tension, not a pressure. 
It is easy to construct flows of arbitrarily low 
energy density which can pass through a mirror 
region, and there is reason to believe that an 


energetic supersonic jet of this type can be stable.°® 


Theoretical analysis of this entire problem 
would be a huge task. Fortunately, certain sim- 
ple experimental expedients should guarantee a 
large degree of trapping. For example, the 
cusped configuration can be made slightly asym - 
metric by misaligning the two coils. The in- 
coming jet will then miss the opposite cusp and 


be trapped. Or else, two streams can be in- 
jected from opposite ends.’ It should be re- 
marked that a certain amount of trapped trans- 
verse magnetic field in the jets is helpful to 
make them collide instead of passing through one 
another as has been observed with “plasmoids.” 
Otherwise one might be able to depend on the 
“two-stream” instability to break up the jets. 
Again, too much symmetry is undesirable since 
it might lead to a spurt sideways through the 
cusped diametral plane, but a slight misalign- 
ment (either in direction or in timing) would 
probably mitigate this effect. 

Generally speaking, trapping is apt to be quite 
efficient in any apparatus with a large approxi- 
mately field-free region (even in a mirror 
machine’), especially if there is no high degree 
of symmetry. It is possible that both the high-, 
and low-8 injection techniques could prove to be 
useful, the first to create a plasma and the 
second to maintain it against cusp losses for the 
required time. 





"The work presented in this paper is supported by the 
U. S. Atomic Energy Commission Computing and Ap- 
plied Mathematics Center, Institute of Mathematical 
Sciences, New York University. 

1J, L. Tuck, Phys. Rev. Letters 3, 313 (1959), 
where this is referred to as “entropy trapping.” 

7H. Grad, Atomic Energy Commission Report TID- 
7520, 1956 (unpublished), p. 148. 

3A concise summary is given by J. Berkowitz et al., 
Proceedings of the Second United Nations International 
Conference on the Peaceful Uses of Atomic Energy, 
Geneva, 1958 (United Nations, Geneva, 1958), Vol. 31. 


‘Hartland Snyder (private communication, 1955). 

‘For example, see H. Grad, Atomic Energy Com- 
mission Report NYO-7969, 1957 (unpublished). 

6D. Finkelstein, G. A. Sawyer, and T. F. Stratton, 
Phys. Fluids 1, 188 (1958). 

"Reference 3, Appendix. The formula quoted there 
is in error. 

8H. Grad, Atomic Energy Commission Report TID- 
7503, 1955 (unpublished), p. 319. See also reference 
3 








°H. Grad, Bull. Am. Phys. Soc. 3, 288 (1958). This 
theory can easily be generalized to include plasma and 
field intermingled. 

‘This proposal, with the two “guns” an integral part 
of the apparatus rather than external, is being carried 
out by a group at Stevens Institute. 
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STUDY OF PHYSICS AND CHEMISTRY OF SURFACES FROM FRUSTRATED 
TOTAL INTERNAL REFLECTIONS 


N. J. Harrick 
Philips Laboratories, Irvington-on-Hudson, New York 
(Received February 3, 1960) 


It has been suggested’ that the physics and 
chemistry of surfaces of optically transparent 
materials, particularly semiconductors, may be 
studied through the analysis of the spectrum of 
totally internally reflected radiation. This seems 
feasible since it is known that the totally reflected 
radiation actually penetrates the surface into the 
rarer medium to a depth given by 


«= 0.693, 
v2 2n(sin6 - n,,?)”*" 





Here \, is the radiation wavelength in the denser 
material, @ the angle of incidence, and »,,, which 
is less than unity, is the relative index of refrac- 
tion.” The penetration depth is about one-tenth of 
a wavelength except near the critical angle; 
where sin? ~n,,, the penetration depth becomes 
large. Goos and Hanchen* demonstrated this ef- 
fect in a remarkable experiment where they ac- 
tually measured a displacement, proportional to 
the penetration depth, of the internally reflected 
beam compared to metallic reflection as shown 
in the insert of Fig. 1(a). 

It is the purpose of this note to give experi- 
mental evidence indicating that infrared radiation 
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FIG. 1. Sample shapes for total internal reflection 
measurements. 
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which penetrates the surface interacts with im- 
purities on the surface at molecular resonance 
frequencies and that this interaction is compar- 
able to what would be observed in a transmission 
measurement for impurity films on the surface 
thinner than a penetration depth. Since there is 
a prohibitive power loss in transmission through 
several separated layers whereas internal re- 
flections are total, the reflection technique has 
the obvious advantage over transmission meas- 
urements where many encounters with the sur- 
face are necessary in order to gain sensitivity. 

The sample‘ used is shown schematically in 
Fig. 1(a) where the angle of incidence is 45 de- 
grees. It was cut, ground, and diamond-polished 
in one piece from 30 ohm-cm, p-type germanium. 
The entrance face had dimensions of about 1/4 
in. X3/4 in. It should be noted that after eight 
reflections the exit beam is axial with the en- 
trance beam, which simplifies the optics of the 
infrared spectrometer. With narrow spectrom- 
eter slits, thin samples can be used and the nun- 
ber of reflections can thereby be increased for 
the same length of sample. Figure 1(b) shows 
how the number of reflections can be increased 
by curving the surface in an appropriate way. 
On each reflection the angle of incidence is in- 
creased. Care has to be taken so that the angle 
of incidence is not decreased beyond the critical 
angle when the beams would be lost. To see if 
this technique might work, known impurities 
were placed on the surface and their spectra 
were analyzed with the aid of a Perkin-Elmer 
Model 12C spectrometer. 

The sample shown in Fig. 1(a) was mounted in 
a cell so that all of the surfaces except the en- 
trance and exit faces could be exposed to a liquid 
or a gas. Figure 2 compares the spectrum of 
the sample with no impurities intentionally added 
to the surface to that observed when water sur- 
rounded the sample. The 3- and 6-micron water 
absorption bands are clearly visible in the latter 
case. When the signal strength of the 3- micron 
band was compared to that observed in a trans- 
mission measurement through a thin film of 
water, the depth of penetration was determined 
to be about 0.15 micron, assuming a comparable 
strength of interaction in the transmission and 
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FIG. 2. Spectra of beam IR of Fig. 1(a). The lines 
attributed to CO, indicate its presence in the atmos- 
phere. 


reflection measurements. This depth of pene- 
tration compares favorably with that calculated 
from theory. 

A series of measurements was next made on 
the 3.4-micron polyethylene line. The polyethyl- 
ene was dissolved in xylene to a dilution of 0.04% 
by weight and sprayed on the surface of the sam- 
ple. The signal from transmission and reflection 
could be directly compared by analyzing the 
spectra from Jy and Jp in Fig. 1(a). The results 
of these measurements are shown in Fig. 3. The 
signal in the reflection beam does not increase 
with film thickness beyond a certain film thick- 
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FIG. 3. Comparison of the signal strength for trans- 
mission and reflection beams, I 7 and IR of Fig. 1(a), 
for the polyethylene absorption line at 3.4 microns. 

The spectrometer slit width was 0.5 mm for these 
measurements. 


ness. The transmission signal, however, in- 
creases with thickness. For a film thickness of 
less than about 0.5 micron, the signal from eight 
reflections is greater than that from two trans- 
missions. When we compare two reflections to 
two transmissions we see that the signal strengths 
are comparable for a film thickness less than 
about 0.16 micron which is also the thickness be- 
yond which the signal strength in the reflection 
measurements does not increase. This is thus a 
measure of the penetration depth of the radiation 
in the rarer medium. Figure 4 shows the pene- 
tration depth of the polyethylene-germanium 
system, whose indices of refraction are 1.5 and 
4, respectively, as a function of angle of inci- 
dence as computed from the formula in the first 
paragraph. We see that depth determined ex- 
perimentally falls right on the theoretical curve. 
These results indicate that the radiation which 
penetrates into the rarer medium does interact 
with the surface impurities at molecular reso- 
nance and that the interaction strength is com- 
parable to that occurring in transmission meas- 
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FIG. 4. The penetration depth versus angle of in- 
cidence for total internal reflection at the germanium- 
polyethylene interface. The experimental point is ob- 
tained from the experimental data of Fig. 3. 
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urements for films whose thicknesses are equal 
to or less than the penetration depth. Further- 
more, this interaction does not represent just a 
power loss at resonance resulting from an in- 
crease in the refractive index in the polyethyl- 
ene, and thus upsetting the condition of total in- 
ternal reflection at the germanium- polyethylene 
interface. If this were the case, the reflected 


radiation would always penetrate the whole of the 
polyethylene film and any radiation not absorbed 


would be reflected from the polyethylene-air 
interface and the signal strength would always 
increase with film thickness. Near the critical 
angle such an effect could undoubtedly be ob- 
served. 

Our results indicate that through the use of a 
more sensitive spectrometer and more reflec- 
tions this technique should indeed work in the 
study of surface chemistry and we expect that it 
should be possible to detect one to ten atomic 
layers. Furthermore, such a technique should 


————_ 


be valuable in the study of thin films. 





‘Discussion by N. J. Harrick following paper by 
R. P. Eischens at Second Conference on Semiconduc- 
tor Surfaces, Naval Ordnance Laboratories, White 
Oak, Maryland, December, 1959 |J. Phys. Chem. 
Solids (to be published).] It was brought to the au- 
thor’s attention after the work described here was 
completed that Dr. J. Fahrenfort of the Royal Dutch 
Shell Laboratories, Amsterdam, has described what 
appears to be a similar technique to observe the spec- 
tra of organic materials on silver chloride at the 
Fourth International Congress on Molecular Spectros- 
copy, Bologna, September, 1959 (unpublished). 
2Julius A. Stratton, Electromagnetic Theory 
(McGraw-Hill Book Company, Inc., New York, 1941). 
3F. Goos and H. Hanchen, Ann. Phys. 1, 333 (1947). 
‘Semiconductor samples similar to these have been 
used in the study of other properties of semiconductor 
surfaces. N. J. Harrick, J. Phys. Chem. Solids 8, 
106 (1959). N. J. Harrick, J. Opt. Soc. Am. 49, 376 
(1959). 








MECHANISM OF ANTIFERROMAGNETISM IN CHROMIUM 


A. W. Overhauser and A. Arrott 
Scientific Laboratory, Ford Motor Company, Dearborn, Michigan 
(Received January 6, 1960) 


An experiment is proposed which may in prin- 
ciple distinguish between three mechanisms of 
antiferromagnetism, one of which may obtain in 
chromium. Corliss, Hastings, and Weiss’ have 
observed that each magnetic superstructure re- 
flection in the neutron diffraction of a Cr single 
crystal is split into six spots equidistant from 
the reciprocal lattice point and which lie along 
[100] type directions from that point. This re- 
sult can be interpreted in several ways: 

A. Cr has basically a Néel type antiferromag- 
netic structure, but with periodic antiphase 
domains." 

B. The localized spins lie in a (100) plane and 
spiral with a fixed period in the direction per- 
pendicular to the plane.’ 

C. The localized spins are oriented by their 
interaction with a static spin density wave in the 
conduction electron gas.* 

Each antiferromagnetic domain (or group of 
antiphase domains for case A) contributes to 
two of the six magnetic reflection spots asso- 


ciated with a given reciprocal lattice point. That 
all six spots occur [and, for example, have equal 
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intensity at the (1,1, 1) point] indicates that the 
unit vector \ characteristic of a given domain 
(as defined below) is parallel or antiparallel to 
one or the other of the three [100] type axes on 
a statistically equal basis. For the three models 
under consideration, i is: 

A. perpendicular to the antiphase domain 
boundaries; 

B. parallel to the axis of the spiral; 

C. parallel to the wave vector of the spin den- 
sity wave. 

An experiment of potentially great interest it 
to cool a single crystal of Cr through the Néel 
point to low temperature in a large magnetic 
field parallel to a [100] crystal axis, and then to 
study the relative intensities of the six magnetic 
reflection spots. If the applied field is sufficient 
to upset the statistical equality of the three [100] 
axes with regard to 3, the relative intensities of 
the spots will be modified, and should allow one 
to distinguish between mechanisms A, B, and C, 
or at least between C and A or B. Such field- 
cooling effects may reasonably be anticipated. 
Magnetic dipole-dipole interactions contribute 
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an energy dependent on the angle between X and 
the axis s of spin alignment. Usually the most 
favorable orientation is x perpendicular to s, as 
is observed to be the case for Cr.’ The spin 
axis $ will in turn be correlated with a magnetic 
field H if the sample is cooled through the Néel 
point with H applied. (It is helpful to recall that 
magnetocrystalline a anisotropy is zero at the 

Néel point, so that s is free to adjust according 
to the direction of H. ) 

Consider now an antiferromagnet conforming 
to model A or B. It is well known that the mag- 
netic susceptibility of sucha system is a maxi- 
mum when s is perpendicular to H. Consequently 
3 will align parallel to one of the two [100] axes 
perpendicular to H. Subsequently, for case A, 
\will be formed either parallel to H or parallel 
to the other [100] axis perpendicular to H. For 
case B, i will be parallel to H since the spiral 
rotation of S sweeps through the plane perpen- 
dicular to H. For mechanism C, however, the 
susceptibility is a maximum when s is parallel 
oH.’ A brief explanation is as follows. The 
axis of quantization of the static spin density 
wave is free to precess around the applied field 
direction. This precession gives rise to a large 
rotating effective (exchange) field H, (in the 
sense of a magnetic resonance experiment) 
which depolarizes the localized magnetic mo- 
ments (similar to the effect of resonance satura- 
tion). As a result, the maximum susceptibility 
occurs when § is parallel to H. Subsequently, 
magnetic dipole-dipole interactions favor x 
being perpendicular to H. The relative intensity 
ofthe two neutron reflection spots parallel to 
Hfrom, say, the (1,1, 1) reciprocal lattice point 
compared to that of the four spots perpendicular 
toH are: 


No field, 1:1; 
field-cooled, A. 2:1, 
B. 1:0, 
c.. Gz. 


The foregoing intensity ratios would allow one 
distinguish between the three mechanisms 
cmsidered. However, these ratios apply only 
{domain orientation adjustment resulting from 
he field-cooling effect is complete. Very likely, 
tagneto-elastic effects associated with domain 
brmation will permit only partial correlations 
tween X and H. In such an event it would be 
iifficult to distinguish between A and B. How- 





ever, one should still be able to distinguish be- 
tween C and A or B since the intensity ratio 
becomes smaller than unity in the former case 
and larger than unity in the latter. 

Whether or not mechanism C operates in Cr 
is of particular interest because some remark- 
able properties of Cr, in addition to the splitting 
of the neutron reflections, can be explained by 
this model. Shull and Wilkinson® have observed 
that the Néel point of a powdered Cr sample is 
45% higher than that of bulk Cr. This result has 
been confirmed by Corliss, Hastings, and Weiss.’ 
Such an effect can be explained as follows. It 
must be remembered that the wavelength of the 
static spin density wave (as determined by the 
splitting of the neutron reflections) differs by 
about 7% from the lattice constant. Consequently 
a significant fraction of the localized spins in an 
antiferromagnetic domain will lie near nodes of 
the wave and will contribute little to the condensa- 
tion energy of the antiferromagnetic state. If it 
were possible for a domain to terminate spatially 
whenever the nodes of the wave begin to fall near 
lattice sites, and for a new domain to begin 
having a more propitious correlation between 
the phase of its wave and the lattice sites, the 
condensation energy and the Néel temperature 
could be increased by as much as a factor of two 
(neglecting domain wall energies, etc.). Con- 
sequently the observed effect is easily accounted 
for, provided one assumes that mechanical im- 
perfections present in the powdered specimens 
provided a mechanism for nucleating (and stabili- 
zing) appropriately small antiferromagnetic 
domains. This explanation seems all the more 
likely since Hastings® has pointed out that the 
neutron reflections would not be split if the do- 
main size were small enough to explain the 
appreciably higher Néel temperature. Indeed, 
the magnetic reflections from his powdered 
specimens were not split.° The behavior of a 
particular sample may well be a function of its 
metallurgical history. 





'L. M. Corliss, J. M. Hastings, and R. J. Weiss, 
Phys. Rev. Letters 3, 211 (1959). 

*T. A. Kaplan, Phys. Rev. 116, 888 (1959). 

°A. W. Overhauser, Phys. Rev. Letters 3, 414 
(1959); J. Phys. Chem. Solids (to be published). 

4A. W. Overhauser (to be published). 

°C. G. Shull and M. K. Wilkinson, Revs. Modern 
Phys. 25, 100 (1953). 
8J. M. Hastings (private communication). 
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NEW MECHANISM OF ANISOTROPIC SUPEREXCHANGE INTERACTION 


Téru Moriya* 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received January 25, 1960) 


Anisotropic exchange interaction as a combined The theory is an extension of the superexchang 
effect of spin-orbit coupling and exchange inter - theory to include the effect of spin-orbit coupling 
action has been discussed by several authors.’ A new formalism recently developed by Ander- 
These previous studies, however, were re- son® is used throughout the calculation. The one. 
stricted to the so-called pseudo-dipolar coup- electron states localized on the cations are ex- 
ling which is of second order in the spin-orbit pressed by Wannier functions and the one-electry 
coupling. The purpose of the present note is to Hamiltonian consists of the energy of the local- 
provide a mechanism for an interaction of the ized states and the transfer terms among them. 
form: We denote the two degenerate wave functions 

D- (S, x3,] ' (1) localized at the position R as follows: 


-R - 
which is linear with respect to the spin-orbit oF ) ad ™ 


coupling, important in magnetic crystals of low 
symmetry, and generally much larger than the 
pseudo-dipolar coupling treated so far. Such an 
interaction was first suggested from purely 
symmetry grounds by Dzialoshinski.” 


where ¢ and 4 mean the spin almost up and al- 
most down, respectively. Because of the spin- 
orbit coupling the spin-up and spin-down states 
are mixed in these wave functions. The one- 
electron terms of the Hamiltonian are written 
as 





E “22 Ze, (Reg (Ray gf) + 2 [yg BBO y (Ria (R)+a “(Ra (R)} 


2 Pus ¥ ’ ™ ’ ’ 
+C (R Ryo 4 (R jor, (R) a (R dar RD} 


¥ Vas . ’ + _ “ , 
+C, (R'-R)a ,, (Ra (R)+C , “(R'-R)a. , (R Jo. , ()], (2) 


, 


n'y 








where 
* 
b  (R’-R)+C *(R'-R)=[¥ (y-R)H.Y (r-R)dr, 
n'n n'n n'h 1 nt 
* 
b  (R’-R)-C *(R’-R)={¥ (r-R)H.Y (r-R)dr, 
n'n n'n n'y 1 ny 
- * 
CR -R)={¥ (r-RYH,Y (7 -R)dr, 
= 
Cc “(ar -R)=fv (y-R)H.© ,(r-R)dr, (3) 
n'n n't 1 nt 
and 


a» (R) and a,4(R) are creation and annihilation 
p? i operators of electrons in the state represented 
H, =5—+V(r) +=—s3S: [gradv(r) xp], by ¥,4(7 -R), etc. The most important contri- 
1 2m 2m*c? n 
bution to the superexchange arises from the 
configuration mixing of polar states due to the 
transfer terms of (2). The second order per- 













Cc. *(R’-R)=C _, *(R’'-R)+iC . >(R’-R). 
n'n n'n n'n 
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turbation gives not only the ordinary isotropic 
coupling (b? terms) but an anisotropic coupling 
coming from bC terms: C* terms are much 
smaller than bC terms and make a contribution 
to the usual pseudo-dipolar coupling. We shall 
consider as an example the case of one electron 
per atom with the ground orbital state assumed 
to be nondegenerate, i.e., the orbital moment is 
quenched. The second order bC terms are cal- 
culated by using the relations: 


. * 
a. (R)a ,(R) “a (R)a (R) - 2s (R), 


a. (R)a (R)=S (R), 


and 
a 
a, (Ra ,(R)=S (R), 


for the occupied (ground) states. The empty 
states do not take part in the second order 
coupling. The result is as follows: 


E, ,,%=D™-[S(R) xS(R’)], 


R,R (4) 


where 


5m =H [o(R -R’')C(R’ -R) -C(R -R’)b(R’ -R)], 
(5) 
the transfer integrals being those between the 
ground states of the ions at R and R’, and U 
being the energy required to transfer an electron 
from an ion to its nearest neighbor, thus making 
apolar state. The transfer integrals b and C 
are calculated by using a perturbation method 
assuming that the orbital level separation in an 
ion is larger than the spin-orbit coupling. The 
bterms are the ordinary transfer integrals with- 
out spin-orbit coupling and the C terms are ob- 
tained, up to the first order in the spin-orbit 





It is interesting to note that all the terms 
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+¢(R’ -R), 
(6) 


where n,n’ represent the ground orbitals and 
m,m’ the excited orbitals, Tinn'®) the matrix 
element of the orbital angular momentum of the 
ion at R, A the spin-orbit coupling constant, and 


c(R’ -R) 


h * - 
-aale,, (ry -R’)[gradV(r) xP]e (7 i. 


yn (7 -R) being the ground orbital wave function. 

From the third order perturbation, which in- 
cludes transfer terms 6} and C between a singly 
occupied and an empty Kramers’ doublet and the 
intra-atomic exchange interaction, we get the 
same type of coupling as (4). This term is 
J/U (J being the intra-atomic exchange energy) 
times as small as the second order term (4) as 
in the case of the ordinary superexchange. 

The direct exchange coupling also has an 
anisotropic part of the form (4). The coupling 
constant is given by 

i *(a)-1, (@) 
_ sf mn "mn 
ex om En) -€ () 
l R’)-1 ’ 
m'n’ ( ) mnt ) 
_ = R’ 
€ AR) -€ (R") 





dJ(nn'mn’‘) 


* 





J(nn nm (8) 


m’ 


where 


’ = “ bi ’ e ’ ey 
J(nn'mm )=fJe, (r, -R)y,, ("5 R rom ? n"9 R)e My -R dt dr. 









unit cell and a center of inversion is located half- 
way between them, this coupling vanishes. When 
the crystal symmetry is low, however, this 
coupling may be the most important anisotropic 
coupling between the spins. Though quantitative 
calculation of this coupling is very complicated, 


treated above which are linear in the spin-orbit 
coupling have the form (1). This coupling van- 
ishes when the symmetry of the crystal is high. 
for example, when there is only one magnetic 

ion in a unit cell and there is a center of inver- 
sion, or when there are two magnetic ions in a 
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we can estimate the order of magnitude in rela- 
tion to the isotropic superexchange. As is easily 
seen, all the three contributions to D treated 
above are roughly (Ag/g) times the correspond - 
ing contributions to the isotropic superexchange 
interaction.‘ Though there may be some reduc- 
tion of the magnitude, depending on the crystal 
structure, this coupling may well be (g/Ag) 
times as big as the pseudo-dipolar anisotropic 
exchange. Extension of the theory to the case of 
many electrons per ion is generally complicated. 
We may expect the above rough order estimation 
is also valid in this case provided the orbital 
moment is quenched. 

The possible existence of this type of coupling 
was first pointed out by Dzialoshinski? from 
symmetry considerations on some trigonal weak 
ferromagnets, a-Fe,O,, MnCO,, and CoCO,, 
though its origin was not clarified. Our mech- 
anism gives the right orders of magnitude to 


explain the magnetic moments of a-Fe,0,, MnCO,, 


and CrF,. 

This coupling is important not only in weak 
ferromagnetism but also in determining the spin 
arrangement in antiferromagnets of low sym- 
metry. As an example we shali consider 
CuCl,-2H,O, one of the most widely studied anti- 
ferromagnets. Because of the low crystal sym- 
metry (space group is V;”) there is a coupling 
of the type (1) between the corner ion and the 
base-center ion; this type of coupling is absent 
between two ions both at the corner (base-center) 
sites. The magnitude of the coupling estimated 
from 4g is about 10% of the superexchange in- 
teraction and D is parallel to the b axis. There- 
fore, the antiferromagnetic spin arrangement 
in this crystal may be different from the one 
accepted at present, i.e., alternating ferromag- 
netic layers in the ab plane with the moment 
pointing in the a direction. We propose instead 
the arrangement shown in Fig. 1. The c com- 
ponent of a spin may be about 10% of the a com- 
ponent. We expect that this spin arrangement 
may give a more consistent understanding of the 
proton resonance data both above and below the 
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FIG. 1. A proposed spin arrangement in antiferro- 
magnetic CuCl,*2H,O. The spins are in the ac plane, 
The black and the open circles represent the ions at 
the positions 0,0, 0 and 3,4,0, respectively. 


Néel temperature’ and the neutron diffraction 
data® for the position of the protons. Further 
details and related topics will be published later. 

I should like to thank Dr. P. W. Anderson for 
his continuous interest in this work. 





* 

On leave of absence from Tokyo Metropolitan Uni- 
versity, Tokyo, Japan. 

1J, H. Van Vleck, J. phys. radium 12, 262 (1951); 
T. Moriya and K. Yosida, Progr. Theoret. Phys. 
(Kyoto) 9, 663 (1953); T. Nagamiya, K. Yosida, and 
R. Kubo, Advances in Physics, edited by N. F. Mott 
(Taylor and Francis, Ltd., London, 1955), Vol. 4, 

p. 1. 

*I. Dzialoshinski, J. Phys. Chem. Solids 4, 241 
(1958). 

3p. W. Anderson, Phys. Rev. 115, 2 (1959). 

44 part of D(2) given by (5) which include ¢(R’-R) 
should be omitted from this order estimation. This 
term is more difficult to estimate than the others. A 
rough order estimation may be given by (spin-orbit 
coupling energy/ionization energy) times superexchange 
energy. 

5N. J. Poulis and G. E. G. Hardeman, Physica 18, 
201 (1952); J. Itoh, R. Kusaka, Y. Yamagata, R. Kiri- 
yama, and H. Ibamoto, Physica 19, 415 (1953). 

8§. W. Peterson and H. A. Levy, J. Chem. Phys. 
26, 220 (1957). 











>a 


Trro- 
ane, 








VoLuME 4, NUMBER 5 


PHYSICAL REVIEW LETTERS 





Marcu 1, 1960 








EFFECTS OF FAST-NEUTRON IRRADIATION ON NICKEL THIN FILMS 


I. Teodorescu and A. Glodeanu 
Institute for Atomic Physics, Bucharest, Roumania 
(Received February 8, 1960) 


In order to obtain more precise information 
concerning the processes and defects produced 
in solids as a result of fast-particle irradiation 
such as those described by Seitz,’ Seitz and 
Koehler,? and Brinkman,’ we have irradiated thin 
films of nickel with fast neutrons. 

The films were evaporated in vacuum (107° 
mm Hg) on glass and quartz substrates, the 
thickness of the films being 200-800 A; the eva- 
poration was carried out at different tempera- 
tures of the substrates in order to obtain differ- 
ent stages of the crystalline structure. During 
the irradiation, the samples were sealed in 
vacuum (10~* mm Hg) or in oxygen. For each of 
the series of the irradiated samples, one test 


sample was put aside for purposes of comparison. 


The irradiation was made for 3.3810" nvt and 
9.45x10"5 mvt at the temperatures T= 59°C and 
T=45°C, respectively. Before and after the 
irradiation, the saturation magnetization and 
coercive force of the samples were measured. 
The structure of the films was studied with the 
aid of the electron microscope and by electron 
diffraction. 

The technique of obtaining the films and the 
methods of studying them have been described 
previously.* 

As a result of fast-neutron irradiation, the 
thin films of nickel were transformed from face- 
centered cubic lattices with a=3.52 A into close- 
packed hexagonal lattices with a= 2.62 A, 
c=4.36 A. This transformation took place in the 
whole film for samples irradiated to 3.38 x10"" 
mt. For samples irradiated to 9.45x10*° nut 
we obtain a mixture of two phases: face-centered 
cubic lattice and close-packed hexagonal lattice. 
The allotropic transformation took place only for 
samples irradiated in vacuum; the samples ir- 
radiated in oxygen were strongly oxidized and 
the oxidized film thus formed prevented the 
allotropic transformation into a hexagonal lattice. 
Such a phase change requires a volume enlarge- 
ment and this is not possible. 

All the irradiated samples show great changes 
in magnetic properties. The magnetic properties 
disappear for the samples which are completely 
transformed into close-packed hexagonal lattices 
(the magnetization falls from 4nl, = 4000 gauss to 


4nl, <6 gauss, and the coercive force falls from 
H,, = 240 oe to H, =0 oe). 

The samples, which are in part transformed 
into hexagonal lattices, show a drop in magneti- 
zation of about 30-50% and their coercive force 
increases to a certain extent; there are cases in 
which it does not vary at all. The increase in 
the coercive force is due to the large stresses 
which appear in the mixture of the phases. 

The samples irradiated in oxygen exhibited a 
drop in magnetization because the proportion of 
pure nickel phase was small and the oxidation 
correction could not be made in the case of 
strong oxidation. That is the explanation of the 
fact that the calculated magnetization is not the 
observed one. The coercive force also grows 
slightly as a result of the stresses. 

The samples irradiated in oxygen showed 
another phenomenon: those samples, which at 
the beginning had a less ordered structure, pre- 
sented after the irradiation a more ordered one 
and the hysteresis cycle was more rectangular. 
The ordered samples were disordered during 
irradiation and the hysteresis cycle was more 
rounded. 

The results of magnetic measurements, elec- 
tron diffraction, and electron microscopy were 
all in agreement. 

In Fig. 1(a) we give the electron diffraction 
pattern for face-centered cubic nickel (nonirra- 
diated). In Fig. 1(b) we give the electron diffrac- 
tion pattern of nickel after irradiation, showing 
that the transformation into a close-packed 





FIG. 1. 
centered cubic (nonirradiated). (b) Close-packed hex- 
agonal (after irradiation). (c) A mixture of the two 
phases. 


Electron diffraction patterns of Ni. (a) Face- 
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FIG. 2. Electron microscope pictures of the same 
samples. 


hexagonal lattice had taken place. In Fig. l(c) a 
mixture of the two phases of hexagonal and cubic 
lattice is shown. 

In Fig. 2 the results of the electron microscopy 
for the same samples are given. 

The results may be explained if we assume the 
following: (1) Thermal spikes were formed and 
their dimensions were ~10° atoms. (2) In the 
region where a thermal spike was formed, the 
lattice lost its meaning as the metal had a tem- 
perature above its melting point. Let us consider 
this as a disordered state. (3) After the thermal 





spike was created, the metal underwent a pro- 
cess of sudden annealing. 

On the basis of these hypotheses, the phase 
transformation may be explained in the following 
way: During irradiation, thermal spikes are 
created followed by sudden annealings. Since 
the annealing period is too short, it is not pos- 
sible to return to the initial state of order (face- 
centered cubic lattice); instead, there is an in- 
termediate state, i,e., the hexagonal close- 
packed nickel. 

The samples irradiated in oxygen do not show 
the transformation phase because their volume 
has been prevented from varying. 

The authors would like to express their grati- 
tude to Miss C. Motoc and G. Alanz for their 
help during the irradiation. 
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FERROMAGNETIC RELAXATION MECHANISM FOR Mz IN YTTRIUM IRON GARNET* 


M. Sparks and C. KittelT 
Department of Physics, University of California, Berkeley, California 
(Received February 8, 1960) 


Yttrium iron garnet (YIG) is an ideal crystal 
for studies of ferromagnetic resonance and relax- 
ation. The outstanding attributes of YIG are 
(A) cubic structure; (B) excellent insulator; 

(C) all magnetic ions identical and in an S state; 
and (D) no detectable magnetic disorder. Spen- 
cer and LeCraw’ have reported ferromagnetic 
relaxation measurements on highly-purified 

YIG at low temperatures. Their work, together 
with that with Fletcher,’ is the first to distinguish 
in a ferromagnet the relaxation of the magnetic 
moment component parallel to the static mag- 
netic field. Their best experimental estimate of 
1/T, at 9.34 kMc/sec between 2 and 30°K is fitted 
roughly, after correction for impurity relaxation, 
by 


1/T, =1x10°T sec", (1) 
with T in deg K; below 5°K the experimental 
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curve flattens out above the values given by (1). 
We propose here an explanation of a significant 
part of the observed relaxation, giving first a 
brief summary of the likely course of the entire 
relaxation process, as well as one can tell at 
present. 

In the absence of special paramagnetic impuri- 
ties, the linewidth is believed* to arise from the 
mixing of the uniform precession (by static sur- 
face and volume imperfections) with spin wave 
modes degenerate with the uniform precession. 
These degenerate spin waves, which we call $ 
modes, then relax by Raman scattering of thermal 
magnons through the magnetic dipole interaction. 
The rate of this process determines 7,, pro- 
vided that the coupling between the uniform and 
S modes is not a bottleneck; at sufficiently high 
temperatures, the bottleneck will exist if AH is 
constant with temperature while 1/7, increases, 
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then 1/7, vs T will appear to flatten. The output 
magnons of the Raman process will tend by 
further Raman processes and by a four-magnon 
exchange interaction process‘ to maintain ap- 
proximately a spin temperature Tsp: One prob- 
ably may not speak of a spin temperature below 
about 5°K for input power ~1 watt/cm*. Except 
at low temperature, the spin temperature may 
be closely equal to the lattice temperature. After 
a study of the magnon-phonon scattering proc- 
esses treated by Kaganov and Tsukernik,°* we 
conclude that the Cerenkov process is dominant 
and is probably adequate down to 5°K. If trace 
amounts of impurities are present which relax 
the uniform mode (and all other modes) even 
slowly, then substantial additional powers can 
be handled, because of the large number of spin 
wave modes. Thus 1 watt/cm®* can be trans- 
ferred from the spin system to the lattice at 5°K 
with a temperature drop ~0.3°K by an impurity 
mechanism giving every mode a 7, =107° sec. 

A failure of spin-lattice contact will flatten the 





ay = -87L(uM | /vy"Dalk, -f, +k 3) 2k 


By standard methods we find, for the relaxa- 
tion time T, of a spin wave k, making an angle 
§ with the z axis, 


3 
167 M (ke gT Je, F(8) 


1 
Dh? w, , (3) 





where F(@) =1+(11/2)sin?6 - (23 /4)sin*6, and D 
defines the high-frequency dispersion relation 
hw=Dk?; w, is the resonance frequency. The 
angle 6 is determined® for a sphere by 


(hw)? = [hy(H - 41M _/3) +Dk,?] 
x[hy(H -4nM _/3)+Dk,? +hy4nM_sin’6]. 
s 1 s (4) 
In the derivation of (3) it was assumed (A) that 
Dk,’?<hw,; and (B) that (hw,)? «<4(kgT)(Dk,?). 

We take for YIG the values D = 0.55 x10-*8 erg 
cm’ and M, =200. We suppose’ that in the actual 
specimen the most important S modes have k, 
2n/d=3x10° cm, where d=2x107 cm is a 
representative grinding particle diameter; for 
this k,, (A) is comfortably satisfied and (B) is 
satisfied for T>0.3°K. For much larger pits 


measured 1/7, vs T curve at low temperatures. 

We have calculated the relaxation rate of the 
S modes for Raman processes in which, by mag- 
netic dipolar coupling between an S mode and a 
thermal spin wave, the z component of the mag- . 
netic moment changes by +2up. The role of 
magnetic dipole coupling in the related problem 
of relaxation averaged over a thermal distribu- 
tion has been treated by Akhiezer.* A calculation 
similar to that below was made by Kasuya® in 
1954, before the line-broadening process was 
understood, so our present note amounts to a 
numerical refinement, on a firmer physical 
basis, of his result. 

The magnetic field’ associated with a . spin 
wave k of magnetization AM is H = -4n(k- AM)k /k?, 
and the magnetic energy coupling one spin wave 
to the medium gives terms bilinear, trilinear, 
and higher in the spin wave operators a", a. 
The bilinear terms may be diagonalized*; the 
trilinear terms split a spin wave in two and vice 
versa. The trilinear terms may be written, 
taking g=2 and wp=lppl, 


“29. a» ta» 
1 (*,) ~~ _ (2) 





the sharp edges on the surface pits will be of 
major importance, but for small pits and local 
imperfections the details are fortunately not 
important. Equation (3) cannot be applied to the 
larger pits without averaging over the surface 
details. 

We calculate 1/T, =0.7410°T,, sec™ on these 
assumptions, in not bad agreement with (1); here 
Tgp is the spin temperature. For linewidth 
caused entirely by point imperfections, one cal- 
culated 1/T, =0.71x10°T,, sec”. We therefore 
see that M, is expected to relax by magnetic 
dipole interaction with thermal spin waves at a 
rate comparable with observation. This process 
will also limit the propagation of coherent spin 
waves which may be generated in the laboratory. 

We are indebted to R. C. LeCraw and R. C. 
Spencer for discussions of the experiments, and 
to R. Loudon for theoretical conversations. 
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ELECTRON NUCLEAR DOUBLE RESONANCE EXPERIMENTS WITH RUBY* 


R. W. Terhune, J. Lambe,f G. Makhov, and L. G. Cross 
Willow Run Laboratories, The University of Michigan, Ann Arbor, Michigan 
(Received January 14, 1960) 


The polarization of nuclear spins throughout a 
crystal by saturation of an electron spin reso- 
nance of paramagnetic ions present as impurities 
has recently been reported.’»? We have observed 
an apparent inverse of this effect in ruby 
(A1,0,:0.05% Cr***), i.e., when the induced 
nuclear polarization is partially removed by 
saturating an aluminum nuclear spin resonance 
transition, a large decrease in the power ab- 
sorbed by the electron spin resonance of the 
Cr*** ions is observed. We have also observed 
effects associated with inducing transitions be- 
tween the hyperfine levels of the chromium -53 
ions. 

Magnetic field modulation at 5 kc/sec was used 
to observe the edge of the (+ 4= -3) electron 
spin resonance absorption line of ruby with a 
microwave power level about 20 db above satura - 
tion. The spectrum shown in Fig. 1 was then 
obtained by scanning the frequency of a low- 
power rf oscillator connected to a single turn 
of wire around the ruby sample. The five lines 
in the 3-Mc/sec region correspond very well 
with the nuclear magnetic resonance spectrum of 
the host aluminum nuclei in the crystal. The 
same dependence of the splitting between lines 
upon the angle of dc magnetic field and line width 
were observed both here and in the nuclear mag- 
netic resonance spectra. 
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FIG. 1. Spectrum obtained by observing the (+4=-4 
electron spin resonance signal from ruby under satura- 
tion conditions and scanning the frequency of an rf 
generator connected to a coil around the crystal. (7 
= 3285 gauss, 9=0°, »v=9160 kMc/sec, T=4.2°K.) 


The lines near 25 and 75 Mc/sec are associated 
with the I- S splitting of the chromium-53 ions’ 
which have a nuclear spin of 3/2 and a natural 
abundance of 9.5%. The other isotopes of chro- 
mium all have zero nuclear spin. The +1/2 and 
+3/2 electron spin states are very close together 
at our operating point. As a result they inter- 
act strongly as the angle of the dc magnetic field 
is varied, causing large changes in the hyper- 
fine splitting (Fig. 2). Also, the second order 
correction to their hyperfine splitting is appre- 
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0 ! 2 3 4 5 6 7 8 9 10 
ANGLE IN DEGREES BETWEEN C AXIS AND MAGNETIC FIELD 
FIG. 2. Experimental data on angular dependence 
of the hyperfine splittings of the electron spin states 
of the Cr®? ion in ruby. The change in frequency with 
angle results from the mixing of electron spin states. 
(v=9160 kMc/sec, T=4.2°K.) 


ciable, leading to the observed triplets. The 
strength of the lines reduced very rapidly with 
departure from zero degrees. A weak triplet 
associated with S= -3/2 was also observed with 
the lines falling at 71.93, 72.15, and 72.37 Mc/ 
sec at zero degrees. It is interesting to note 
that while all the Cr*®* data were taken by ob- 
serving the S=-1/2 to +1/2 transition, saturation 
of the hyperfine levels in the S=+3/2 state 
caused the largest changes in the electron spin 
resonance signal. 

The strong effect of nuclear polarization on the 
electron spin resonance can easily be observed 
in ruby without applying rf. At microwave power 
level well above saturation, if one rapidly moves 
from one part of the resonance line to another 
a large transient with a decay time of about 5 
seconds is observed. The same decay time is 
observed for the transient following removal of 
tf power sufficient to saturate one of the alumi- 
tum nuclear transitions. Further, the observed 
decay time of the aluminum nuclear polarization 
isalso 5 seconds while 7, for the electron spin 










resonance transition is approximately 0.2 second. 





The decay time of the effects associated with the 
chromium-53 nuclei was also near 5 seconds. 

The effects of nuclear polarization upon the 
electron spin resonance signal were only obser- 
vable when microwave power levels of the order’ 
of or greater than that needed to saturate the 
electron spin transition were used. When the 
electron spin resonance signal was observed at 
low microwave powers, the application of rf had 
no effect upon the signal. Further, when the 
signal was observed during the nuclear polariza- 
tion relaxation time following a quick reduction 
in the microwave power level from above to well 
below saturation, only a 0.2-second time con- 
stant transient was observed. The application of 
rf power had no effect upon the transient, indi- 
cating that 7, is independent of the nuclear 
polarization. 

A decrease in microwave power absorption 
with the application of rf power was always ob- 
served. Asa check, some observations were 
made using amplitude modulation of the micro- 
waves rather than magnetic field modulation. 
Also, the amplitude of the magnetic field modula- 
tion was varied over a wide range with no signifi- 
cant effects being observed. 

It does not seem that the effects reported here 
can be explained by any of the mechanisms pro- 
posed by Feher‘* to explain his observations in 
doped silicon. The slow relaxation after the 
removal of rf power would seem to indicate that 
the host nuclei play an important role in either 
the induced transition probabilities or relaxation 
mechanisms between electron spin states. If 
such were the case, in order to get an appreci- 
able percentage effect in the electron system by 
saturating the nuclei, one would have to reduce 
the relative electronic polarization (N;/N;) to the 
same order as the nuclear polarization (n,/n;), 
i.e., Nj/Nj=np/nj. This idea is consistent with 
our experimental results as the nuclear polari- 
zations are increased by about a factor of’ 40 
and the electronic polarizations are decreased 
by about a factor of 100 under our experimental 
conditions, making the two comparable. 

We have also studied the effect of rf upon 
maser action in ruby,® and have observed a num- 
ber of additional very weak lines. 

We wish to thank Professor C. Kikuchi for 
many helpful discussions. 
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THERMAL EQUILIBRIUM BETWEEN F AND M CENTERS IN POTASSIUM CHLORIDE 


C. Z. van Doorn 
Philips Research Laboratories, N. V. Philips’ Gloeilampenfabrieken, Eindhoven, Netherlands 
(Received February 3, 1960) 


Various models have been proposed for the M 
centers in alkali halides. According to Seitz’ an 
M center consists of an F center associated with 
a neutral vacancy pair. A modification, having 
inversion symmetry, was put forward by Knox.” 
Van Doorn and Haven® have proposed a model 
consisting of a pair of associated F centers (F, 
model). From these models it is clear that the 
formation of M centers from F centers requires 
one F center for each M center formed in the 
case of the Seitz-Knox model and two in the case 
of the F, model. 

It might be possible to distinguish between 
these cases by investigating (a) the kinetics of 
the reaction F~M and (b) the thermal equilibrium 
F=M. For the Seitz- Knox center the equilibrium 
can be written as 


F + vacancy pair=M, 
K =[M]/[F ][vac. pairs]; 


and for the F, center as 
K'=([M)/[FF. 


In the first case the M concentration will be pro- 
portional to the F concentration at a fixed tem- 
perature (the concentration of vacancy pairs 
being constant), whereas in the second case the 
M concentration will vary quadratically with the 
F concentration. 

Theisen and Scott* found the thermal bleaching 
of F centers at about 300°C to be a second-order 
reaction, indicating that two F centers combine 
to form a new center. The thermal equilibrium 
between F and M centers was studied by Smith 
and Scott.° They colored a KCl crystal inhomo- 
geneously by heating for a short time in potas- 
sium vapor in order to obtain a range of F con- 
centrations. It was assumed that the attainment 


2F=M, 
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of the equilibrium between F and M centers is 
rapid compared with the diffusion of the centers 
in the crystal. The F and M concentrations in 
different regions of the crystal were measured 
optically, using a thin beam of monochromatic 
light. A linear dependence between F and M 
centers was found up to a concentration of 10" F 
centers/cm*. This experiment is, however, not 
conclusive because (a) no true equilibrium ex- 
isted and (b) the high M/F ratio found indicates 
that the crystals were not properly quenched. 
Therefore, the following experiment was per- 
formed. KCl was purified by crystallizing it 
twice according to the Kyropoulos method. From 
the single crystal thus obtained, plates were cut 
having the dimensions 4x5x1.3 mm*. These 
were colored by heating to 697°C in potassium 
vapor of various pressures in a nickei tube. A 
heating time of 10 minutes was found to be suffi- 
cient to obtain equilibrium, as no further in- 
crease in optical absorption was observed with 
longer heating times. After quenching in carbon 
tetrachloride to room temperature, which was 
done in the dark to avoid photochemical reac- 
tions, the relative F and M concentrations were 
determined by measuring the extinction of the 
absorption bands at 77°K. With the apparatus 
used for coloring (to be described in a forth- 
coming publication) it took less than one second 
to transfer the crystal from the furnace to the 
quenching bath, thus ensuring freezing- in of the 
high-temperature equilibrium. All crystals had 
the same dimensions in order to obtain equal 
quenching rates. 

The results are shown in Fig. 1. It is seen 
that the M concentration varies quadratically 
with the F concentration, which is inconsistent 
with the Seitz-Knox model but agrees well with 
the F, picture. 
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FIG. 1. Double-logarithmic plot of the concentra- 
tions of F and M centers in equilibrium at 697°C, ex- 
pressed as absorption constants at 77°K of the F and 
M bands, respectively. The number of F centers per 
em? is calculated from the Smakula formula with 
f=0.81 and a half-value width 0.195 ev. 


A rough determination of the equilibrium con- 
stant K’ at other temperatures showed the tem- 
perature dependence to be small, corresponding 
to a heat of formation of the M center of about 
0.01 ev. 

I would like to thank Mr. A. T. Vink for his 
assistance. 
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NUCLEAR RESONANCE FLUORESCENCE IN Au?®’t 


D. Nagle, P. P. Craig, J. G. Dash, and R. R. Reiswig 
Los Alamos Scientific Laboratory, University of California, Los Alamos, New Mexico 
(Received February 1, 1960) 


We have observed the Mossbauer effect’ in the 
77-kev gamma-ray line of Au’®’. A resonance 
curve from which the level lifetime may be de- 
duced was obtained by moving a Pt*®” source 
relative to a Au absorber. Resonance was found 
using as sources either the 8~ parent, Pt’®’, or 
the electron capture parent, Hg’®’. We also find 
that the effect can be either enhanced or dimin- 
ished by placing the source and/or absorber 
atoms in suitable crystal lattices. All measure- 
ments were performed at a temperature of 4°K. 

The shape of the resonance line was determined 
using a Pt*®? source prepared by neutron activa- 
tion of 60 mg of natural Pt. Pulse-height spectra 
taken at intervals for several days following 
activation showed a prominent gamma line at 
71 kev, decaying with the mean life of about 26 
hours, which is characteristic of Pt’*’. Figure 
1 shows the counting rate (corrected for source 
decay during measurement) versus the velocity 
of the source relative to a 200-mg/cm? Au ab- 
sorber. No distinction between positive and 
hegative velocities was made. The velocity half- 
Width obtained is 3.54 0.5 mm/sec. No evidence 
of subsidiary resonances was seen up to the 
highest velocity used, 24 mm/sec. The total 
resonance absorption effect is 5.54 0.5% after 
@ correction (not shown in the figure) is made 
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FIG. 1. Nuclear resonance absorption vs the mag- 
nitude of the source velocity in the 77-kev gamma ray 
of Au!®?, 


for background in the gamma-ray spectrum. 
Using a source prepared by activation of en- 
riched Pt (65.9% Pt'**), measurements of reso- 
nance absorption were made on three Au ab- 
sorber foils of thickness 100, 200, and 400 mg/ 
cm’. The net absorptions were, after correction 
for background (including the unresolved 60-kev 
x ray), 4.020.6%, 7.7+0.5%, and 11.5+0.7%, 
respectively. The velocity curve was similar 
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to that of Fig. 1. Ascribing this nonexponential 
dependence upon absorber thickness to prefer - 
ential absorption of the center of the line, we 
determine the probability f for emission of 
photons in no-phonon processes, and the prob- 
ability f’ of resonant absorption from an incident 
beam of uniform energy spectrum,’ namely, 

f =2347% and f’=3+41%. The calculation of f’ 
is based upon a resonance cross section which 
includes the value*® 3.5 for the internal conversion 
coefficient of the 77-kev Au’®’ line. An earlier 
value* of 5.9 would raise f’ to 4.64+1.5%. 

f and f’ are intimately related to the phonon 
spectra of the source and absorber. In the limit 
of low temperatures, Visscher® has shown that 
for monatomic isotropic crystals this relation 
takes the simple form of a weighted integral over 
the phonon spectrum. We feel that values of these 
weighted integrals are the major results of trans- 
mission measurements using the Mossbauer 
effect. We define an effective temperature, 


o- [eta fmaeroi] 


where g(v)dv is the number of phonon states in a 
frequency range dv, N is the number of atoms in 
the solid, and k is the Boltzmann constant. Then 
for the source, @ is related to f by @=R/in(1/f), 
where R is the classical recoil energy of a source 
atom of energy E and mass m, Rk =E*/2kmc’. 

The same relation relates the absorber phonon 
spectrum to /’. 

We find that for Au, @=54+t§ °K; for Pt, @ 
= 127*27 °K. Our earlier results* on Os and Ir 
yield @=(200+17)°K for Os and @=(1744 20)°K 
for Ir. 

The simplified Debye model of phonon excita- 
tions in solids yields @= 26, where @ is the 
Debye characteristic temperature. On this model, 
@=123°K for Au and @=155°K for Pt, using @ 
=185°K for Au and @=233°K for Pt,° in significant 
disagreement with our results. However, one 
should not place too much reliance on the pre- 
dictions of the Debye model in this instance. 

By using the above transmission results, the 
resonance curve of Fig. 1 may be corrected for 
broadening caused by selective absorption in the 
200-mg/cm? absorber. For this absorber the 
corrected mean life may be written’ 


rat 1238 ’ Sx2.s1(*250) musec, 


where E is the gamma-ray energy (77 kev) and 
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u is the velocity half-width in mm/sec. Using 
our determination of f’ we find 7 =0.82+0.15 
mpsec. The mean life obtained by conventional 
techniques is 1.9+0.2 musec.* 

There is some indication of hyperfine struc - 
ture which is less than the natural linewidth and 
which is slightly asymmetric about zero velocity. 
We are at present investigating this point. Such 
an asymmetry would broaden the observed reso- 
nance curve and decrease /’ below the value 
expected from the Debye model, in agreement 
with our observations. 

In an effort to understand the influence of the 
lattice binding upon the resonance, several ex- 
periments were performed in which the source 
or the absorber was modified. An absorber was 
prepared of the compound AuAl,, containing 
208.5 mg/cm? of Au. The Lindemann formula® 
was used to compute roughly the Debye tempera- 
ture of this compound. With the enriched Pt 
source, a nuclear resonance effect of about 21% 
was to be expected. We were able to observe 
only a marginal effect of about 0.3 + 0.3 %. 

The source lattice may be drastically altered 
through the use of the electron capture parent, 
Hg'*’. Such a source was prepared by neutron 
activation of natural Hg. On the Debye model, 
taking 6 =69°K, the total resonance absorption 
expected with a 200-mg/cm? Au absorber is 
1.1%. The spectrum has, however, unresolved 
60-kev x rays from K capture in Hg which we 
estimate to be 2.5 times as intense as the 77- 
kev line of interest. The effect observed was 
(0.14+0.12)% before correcting for x rays, and 
after correction (0.49 + 0.42) %. 

The small effect in this case may be attributed 
to a low value of the effective temperature ¢@ 
for the mercury, and to a lack of symmetry be- 
tween the crystallographic environments of 
emitting and absorbing Au’®’ nuclei. According- 
ly, we tried to enhance the absorption by adding 
active mercury to compacted gold powder in the 
proportion 10% Hg and 90% Au. The amalgam 
was then diffusion annealed for 70 min at 430°C 
under a helium atmosphere. One day later the 
resonance effect was measured and found not to 
differ from the result with pure Hg. Previous 
X-ray work had shown that samples of this type 
were homogeneous after 24 hours at 430°C. 
Therefore the sample was given an additional 
24 hours heat treatment at 430°C. The resonance 
was now found to be (0.464 0.04) % before back- 
ground correction and (1.6+0.5)% after correc- 
tion. It appears that placing the emitting atoms 
in a lattice which is mostly Au has resulted ina 
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substantial enhancement of the resonance. 

We wish to thank A. D. McGuire for help in the 
early portion of this experiment. We would like 
to thank M. Goldhaber for a helpful discussion. 
We take pleasure in acknowledging the help of 
the Omega reactor group in preparing sources. 
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Van Hove’ has shown that the partial differen- 
tial cross section for the coherent scattering of 
neutrons by an assembly of atoms is related to 
the space-time correlation function, G(r, #), by 


&o a® k’ 
Wee Ahk Er&), (1) 
I'(k, w) = ff expli(k-r - wt)]G(r,t)drdt, (2) 


where a is the appropriate scattering length, k 
and k’ the wave vectors of the incident and out- 
going neutron, hx the momentum transfer, and 
€=hw the energy lost by the neutron. The spa- 
tial part of the transform (2), y(k,#), is given by 


N 
| y(k,t)=N" 2 _ exp[-ik-r(0)] 
1,j=1 t 


x explix-F,(4)]) 7, (3) 


where (-++)~ denotes a thermal average over the 
states of the scattering system, and F(t) is the 
Heisenberg operator for the position of the jth 
atom at time ?. 

Now it can easily be verified, by writing the 
thermal average and r;(¢) in terms of the Hamil- 
tonian of the system, that y(k, t+ ih/2k pT) is a 
real function of ¢. Let this function be denoted 
by y’(k,#). The Fourier transform of y’(k, t) 
with respect to ¢ then defines a function I'’(k, w) 
which is related to I'(k, w) by 


I'(k, w) =P(k, w) exp[hw/2k 7]. (4) 





From the properties of y’(x, ¢) it follows that 





SPACE-TIME CORRELATION FUNCTION FORMALISM FOR SLOW NEUTRON SCATTERING 


P. Schofield 
Theoretical Physics Division, Atomic Energy Research Establishment, Harwell, England 
(Received January 29, 1960) 


























I''(k,w) is even in both k andw. Thus (4) ex- 
presses the condition of detailed balance, which 
we know a priori must be satisfied by the scat- 
tering cross section. 

By inverting the transform (2) and using (4), it 
is possible to obtain a dispersion relation re- 
lating the real and imaginary parts of G(r, ?). 
This may be written in integral form: 


9{G(r, t)}=(k pr/h) 


xf csc(nk Tt’ /h)a{G(r, t-t’)}at’, (5) 


- co 


or, more conveniently, in operator form: 
s{c(, n}=-tan . aye, )}. 6) 


This can be regarded as a constraint imposed 
upon the correlation function as a consequence 
of the thermal equilibrium of the system. 

Vineyard? has suggested that in many systems 
it is possible to replace G(r, ¢) in (2) by its clas- 
sical equivalent G,(f,t). Since the classical 
correlation function is real, the resulting cross 
section does not satisfy the detailed balance con- 
dition. The above analysis suggests that G(r, t) 
should be used to define ['’(k,w), rather than 
I'(k,w), so that this condition is satisfied. Al- 
though in the limit A-0, »’(k, t) and y(k, t) are 
identical, since the variable of transformation 
is the energy transfer divided by #, the trans- 
forms are not identical. 

The use of the classical correlation function in 
this way takes some account of the recoil of the 






239 






















VoLuME 4, NUMBER 5 


PHYSICAL REVIEW LETTERS 











Marcu 1, 1960 








scattering atom, and, as may be seen by con- 
sidering the simple examples of the monatomic 
gas and the isotropic oscillator, gives a cross 
section which is formally correct to first order 
in h. 

Another case discussed by Vineyard? is that in 
which the atoms diffuse according to the Langevin 
equation. By means of the real function y’(x, ¢) it 
is possible to obtain an explicit expression for 
the cross section: 


Wey 2 aS (=) mel 
P(k,w)= > ¢ pas —— a 
a eS 1 
x [Ca+n)*+ 8° (a+n+1)? + B? ’ (7) 








where 7 is the coefficient of viscous damping, 


and 
a =k ,TK?/Mi, B=w/n. 


The result obtained if the simple diffusion 


equation is used to derive G(r, t) corresponds to 
the first part of the first term in (7) so that un- 
less both a and g are small, this approximation 
is not valid. In fact, (7) gives a width for the 
scattered energy distribution less than that given 
by both the simple diffusion model and the gas 
model. This is in agreement with the observa- 
tions of Pelah et al.* on liquid lead. However, 
(7) does not give a width sufficiently small to ex- 
plain the results of Turberfield* on lead with cold 
neutrons. These experiments indicate that the 
atoms of the liquid must make many vibrational 
oscillations between diffusive jumps, and the 
results cannot be interpreted on a basis of con- 
tinuous diffusion. 
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LEAD K ABSORPTION EDGE FOR »-MESON MASS DETERMINATION 


Alan J. Bearden* 
Department of Physics, The Johns Hopkins University, Baltimore, Maryland 
(Received January 26, 1960) 


The recent accurate determination of the 
fulf p ratio for the 1 meson’ and the compari- 
son of the magnetic moment of the 1 meson 
based on this measurement with the value pre- 
dicted by quantum electrodynamics’ depends 
critically on the mass of the » meson. It has 
been shown experimentally that the 3D-2P tran- 
sition in p -mesonic phosphorus lies near the 
Pb K edge energy.® In fact, when considered in 
detail, it appears that the fine structure split- 
ting of the 3D-2P transition in phosphorus will 
cause the lines to straddle the Pb K edge energy. 
Therefore, it seems possible that one may use 
the variation of the mass-attenuation coefficient 
in lead as a function of energy in the K edge 
region as a very accurate measurement of the 
mesonic x-ray energy. The detailed shape of 
the lead K edge must be known precisely both 
in energy scale and mass-attenuation coefficient. 
Knowing these transition energies, and account- 
ing for corrections due to vacuum polarization, 
nuclear size effects, etc., it would be possible 
to make a precise determination of the 1 -meson 
mass. The transition energies for the mesonic 
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phosphorus lines under discussion are given by* 


3Dyj2-2Py2: 425.42 m ,/m,, relative intensity =1; 
3D yz -2Py2: 425.65 m ,/m, relative intensity =9; 


3Dy2-2Py2: 427.75 m /m, relative intensity =5. 


From the quoted “best value” of the 1 -meson 
mass,” m,,/mg = 206.86+0.11, the 3D, -2Py, 
transition would lie at 88041 ev. The best esti- 
mate of the inflection point of the Pb K edge is 
at 88 014 ev,°® and its width is of the order of 100 
ev. The intent of this Letter is to present then 
an accurate measurement of the Pb K absorption 
edge. 

A two-crystal spectrometer using the 2023 
planes of quartz as diffracting surfaces was con- 
structed. Angles were measured by use of the 
precision circle of the Societé Genevoise goniom- 
eter of the Johns Hopkins x-ray laboratory.’ In 
order to gain comparative freedom from scat- 
tering, distances of one meter were used be- 
tween x-ray source, first crystal, second crys- 
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tal, and detector. With this arrangement the 
signal to background ratio was never worse than 
one hundred to one. A commercial high kilo- 
voltage x-ray tube* was used as an x-ray source, 
and a Nal(T1) scintillation counter with single- 
channel pulse-height discrimination detected the 
radiation. The full width of the (1-1) rocking 
curve at the Pb K edge wavelength was 5.5’’. 

The use of this narrow window allowed the de- 
termination of the edge shape without any appre- 
ciable interference due to window width. The 
“thickness effect” as previously reported®’?® is 
negligible in this case, as was shown by meas- 
uring the mass-attenuation coefficient near the 
high side of the edge and varying the thickness 

of the absorbing foil. No change of mass-attenua- 
tion coefficient of more than two percent was 
detected over an absorber thickness from one- 
quarter to five absorption lengths. This absence 
of “thickness effect” is a distinct advantage of 
the high-resolution two-crystal instruments, 

and is due to the smaller unbent crystal diffrac- 
tion widths and the larger dispersion obtainable 
by the two-crystal method. 

The crystal spacing was determined by meas- 
uring the diffraction angles for Mo Kq1, Ag Kg}, 
and W Kg j radiation. The lattice spacing for 
the 2023 planes in quartz at 25°C was found to be 
(1372.117+ 0.010) xu. 

The Pb K edge determination was made with 
the diffracted x-ray beam passing through the 
crystals. As a check, the W Kq line was meas- 
ured with: (1) the diffracted beam passing 
through the crystal; and (2) the diffracted beam 
in “reflection.” No significant difference in dif- 
fraction angle could be detected, once the dif- 
fraction angle in case (2) was corrected for re- 
fraction in the crystal. 

The measured shape of the Pb K edge for 98% 
purity lead foil’* appears in Fig. 1. The energy 
scale is known to +3 ev and the value of the mass- 
attenuation coefficient at any point is stated with 
a probable error of two percent. Since the 
3Dy, -2P,. transition probably lies entirely on 
the high absorption side of the Pb K edge, the 
mass-attenuation coefficient of lead was meas- 
ured at 88 450 ev and found to be (7.30+ 0.1) cm?/ 
gram. 

As the two-crystal spectrometer measures 
the mass -attenuation coefficient for an extremely 
“good geometry,” the mass-attenuation coeffi- 
cient includes both photoelectric absorption and 
scattering cross sections. For any particular 
mesonic x-ray experiment, it is necessary to 








1 


s 
T 


7 
r Me ten? sqrem) 


w 
ye 


nN 
T 








87 650 87 900 87 350 88 000 68 050 88 100 68150 


Electron Volts 


FIG. 1. PbK absorption edge. 


ascertain the amount of scattered radiation which 
enters the detector in order to arrive at the ap- 
propriate mass-attenuation coefficients for the 
unavoidable “poor geometry.” Unfortunately, 
present measurements of the attenuation of 

3D -2P mesonic x rays of phosphorus in lead® are 
not suitable for use in a precise determination 

of the ,»-meson mass. Therefore, new experi- 
ments are essential to determine with high preci- 
sion the attenuation in lead for these mesonic 
X-ray transitions. 

The author is indebted to Dr. R. L. Garwin and 
Dr. S. Penman for discussions of this problem, 
and to the Nevis Cyclotron Laboratory of Colum- 
bia University for supplying the high-voltage x- 
ray tube. Many helpful discussions as to the 
measurements were had with Professor L. Madan- 
sky and Professor J. A. Bearden. 
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1* -p TOTAL CROSS SECTIONS IN THE RANGE 450 Mev TO 1650 Mev* 


Thomas J. Devlin, Barry C. Barish, Wilmot N. Hess, Victor Perez-Mendez, and Julius Solomon 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received February 8, 1960) 


A series of counter experiments is being car- 
ried out at the Berkeley Bevatron and 184-inch 
cyclotron to measure total and differential scat- 
tering cross sections for pions at energies above 
the first resonance.’ In this Letter we report 
the results of measurements of the total cross 
section for positive and negative pions on pro- 
tons in the energy range from 450-Mev to 1650- 
Mev pion kinetic energy (lab). 

The experimental arrangement is shown in 
Fig. 1. The pions were produced in an internal 
Pt target, plunged into the circulating proton 
beam of the Bevatron. The Pt target was located 
in the field-free tangent tank, allowing us to 
examine both positive and negative pions with the 
same experimental setup. Our magnet system 
consisted of an 8-in.-bore focusing quadrupole 
doublet, an 18 x36-in. bending magnet which 
deflected the particles through an angle of 23 deg, 
and a 4-in.-bore focusing quadrupole triplet, 
used as a field lens. The absorber was a 48-in. 
target of liquid hydrogen.* 

The beam incident on the liquid hydrogen tar - 
get was monitored by coincidence of a set of 
plastic scintillation counters, M (1.5 and 2 in. 
in diameter), and a gas Cerenkov counter,° for 
discrimination against protons in the positive 
beam. In the energy range in which protons of 
the same momentum had a velocity less than 
0.8c, a liquid nitrogen Cerenkov counter was 
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used in place of the gas counter. The trans- 
mission of pions was measured by a set of scin- 
tillation counters, S, to S,, in coincidence with 
the monitor telescope. The scintillation counter, 
S,, was placed in coincidence with S, to S, to 
reduce the accidental background. Counters S, 
to S, subtended different sclid angles from the 
target, ranging from 1.5 to 7.1 millisteradians. 
The cross section was determined for each of 
these counters and extrapolated to zero solid 
angle. Eight transmission measurements were 
made at each energy, four with the hydrogen 
target empty and four with it full. 

The momentum of the pion beam was deter- 
mined by deflection through the bending magnet. 
It was measured by the current -carrying-wire 
method at four different times, before, during, 
and after the experiment. The consistency of 
these measurements was of the order of 1%. 
The momentum was also measured with a gas 
Cerenkov counter pressure curve. This gave 
the same result, but with somewhat less ac- 
curacy. The momentum spread determined by 
the counter telescope and the uncertainty in 
momentum were both 1.5%. 

It was necessary to make several corrections 
to the data. Several species of accidental counts 
were monitored continuously during the experi- 
ment. These were negligible for negative pions. 
A fairly significant accidental rate was observed 
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in the monitor telescope for positive pions. This 
was due to the high flux of protons in the beam, 
and can be corrected for with the known proton- 
proton cross section. 

The beam is also contaminated by muons from 
rdecay, and by electrons. The muon contamina - 
tion originates from two regions: before the 
pending magnet, and after it. The contribution 
from pion decays before the bending magnet was 
measured directly by means of a gas Cerenkov 
counter pressure curve. An exact calculation 
of this contribution is difficult, and we were 
forced to make some simplifying approximations. 
The measured contamination agreed very well 
with the approximate calculations. The contri- 
bution from decays after the bending magnet was 
not measurable by the pressure-curve method, 
but it is calculable in a simple straightforward 
manner. The total muon contamination varied 
from 2.54+1.0% at 1650 Mev to 9.5+1.0% at 450 
Mev. The electron contamination was also meas- 
ured by the pressure curves, and it was esti- 
mated to be less than 1%. 

The effect of Coulomb scattering in the hydro- 
gen target was taken into account in the extra- 
polation of the cross section to zero solid angle. 
Athigher energies, the effect was negligible at 
all solid angles; at lower energies, the smaller 
solid angles were affected. Where necessary, 


the cross sections were corrected by a method 
similar to that of Sternheimer.® 

The number of counts taken in each measure- 
ment gave a statistical uncertainty of 1% or less 
in the cross section at each solid angle. The 
fluctuations in the series of eight transmission 
measurements for each particular cross section 
were consistent with this. The errors on the 
points in Fig. 2 include the systematic errors 
due to muon, Coulomb, and accidentals correc- 
tions in addition to counting statistics. 

The total cross sections for positive and nega- 
tive pions on protons are shown in Fig. 2. Our 
results show peaks in the 7~ cross section at 
600 +15 Mev and at 900215 Mev. These are 
somewhat lower energies than indicated earlier 
by the MIT experiment,’ but they are in substantial 
agreement, within statistics, with the measure- 
ments of the group at Saclay,® which are also 
plotted in Fig. 2. The existence of a maximum 
in the 7+ cross section was indicated earlier by 
an experiment at Brookhaven.® Our results, in 
this and in a previously reported experiment,’ 
confirm its existence and show it to be centered 
about 1350 Mev. 

We wish to acknowledge our indebtedness to 
Professor Burton J. Moyer and Professor A. C. 
Helmholz for their constant help and encourage- 
ment, to Michael Longo and John Atkinson for 
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FIG. 2. Total *-p cross sections vs pion kinetic energy. 
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their valuable aid in all phases of the experiment, 
and to the Bevatron staff and crew. 
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QUESTION OF THE EXISTENCE OF A STRANGENESS 2 MESON* 





D. J. Prowse 
Department of Physics, University of California, Los Angeles, California 
(Received February 3, 1960) 


Recently Yamanouchi’ has suggested that there 
is strong evidence for the existence of a strange- 
ness 2 meson with a mass of about 720 Mev. 
Among several events enumerated by Yamanouchi 
in support of this contention is the Bristol anom- 
alous K*-meson decay reported by Prowse and 
Evans.” 

Yamanouchi points out that if the identity of the 
primary particle in this event is unknown one 
could set up a scheme, D+~7++K° rather than 
Kt ~a*++7°+y or Kt —2*+x°, because the identi- 
ties of all the neutral particles involved are un- 
known. The argument on whether or not the 
decay involves two or three bodies is based upon 
the fact that in a similar event found by the Colum- 
bia group,® the 7+-meson energy was the same 
(approximately 60 Mev). If the existence of two 
decays giving the same 7*-meson energy can be 
considered as possible evidence of a unique 7* 
energy, then one is justified in postulating the 
decay mode K*~2*+x° or, if the nature of the 
primary particle is undetermined, D* -1*+K° 
also. The purpose of this Letter is to point out 
that the mass of the primary particle is known 
in the Bristol event and it is not consistent with 
the mass value of 720 Mev. 

The event in question was found during a system- 
atic study of the interactions of K* mesons with 
nuclei.* In the course of this work it appeared of 
interest to determine the decay modes of those 
K* mesons which had interacted and to compare 
them with the decay modes of those K* mesons 
which had not interacted.* The decay particles 
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of the K* mesons which had been followed to rest 
in the search for interactions were therefore 
examined. Details of the pickup criteria used 
for finding and following these tracks are given 
in reference 4; the details of the measurements 
on the secondary particles are given in reference 
5. In the course of the grain density measure- 
ments on these secondaries it became important 
to normalize the values obtained very carefully 
in each pellicle utilized. The normalization use 
was the grain density of the beam 7 mesons 
(1.01 Xminimum); this value was carefully de- 
termined in each plate as outlined in reference 5 
and shown in Fig. 2 of reference 2. A number of 
K*+-meson identities were established by g” vs 
measurements because not all secondaries were 
visible owing to the low minimum grain density 
in the stack. In particular the identity of the 
primary of the anomalous event was established 
in this way. Grain density measurements were 
made at three residual ranges and the values 
were normalized using the minimum -ionization 
results mentioned above. The K*-meson track 
was quite flat (1.4 cm/pellicle) and so no diffi- 
culties were encountered because of steepness. 
The three values obtained are shown graphically 
in Fig. 1. The expected variations of g* vs R 
are shown for 7 mesons, K mesons, D particles, 
and protons. The mass of the primary is thus 
determined to be 532+24 Mev. The identity can 
be seen to be well established. 

The high strangeness number of 2 suggested by 
Yamanouchi for the new particle is sufficient to 
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FIG. 1. Grain density measurements plotted against 
sidual range (R). The expected variations of g* 
normalized to minimum) with R for the various parti- 
es are computed from the curves of G. Alexander 
 R.H.W. Johnston [Nuovo cimento 5, 363 (1957)] 
org" values below 2.9. Above this value the curves 
e computed from the values of (e")p. p, given by 
?.H, Fowler and D.H. Perkins [Phil. Mag. 46, 587 
1955)]. 


xplain the apparent absence of conclusive proof 
{its existence. It is likely to be removed from 
ny mass-separated beam which has been used 
eretofore and the energies of the separated 
trange-particle beams have not been high enough 
»produce it. The only experiments which 

night have detected the particle are those in- 
olving momentum-separated beams only. The 
mulsion stack in which the Bristol anomalous 
\*meson was found was exposed to such a beam 
oi positive sign) of 375 Mev/c momentum. Dur- 
ug the systematic scan for K* particles, a 

search was also made in this stack for stopping 
articles of mass 500m, and mass 1480mg; this 
search was inspired by previous reports of parti- 
les of this mass. No particles were found stop- 
ing and decaying in the regions of the emulsion 
lack corresponding to these mass values. If 

he particles have the same lifetime as the K* 
meson, the ratio of (D+ or X*)/K* at production 
sless than 0.5%. 

The evidence for the existence of the particle, 
*, is partially listed by Yamanouchi. The re- 
bulls given in this paper show, however, that 
he 2 anomalous events® have to be deleted from 
his list. This is unfortunate for two reasons: 
4)The exact mass value is not now known, and 
) it materially decreases the probable lifetime 





of the particle. Both of these factors will make 
it harder to detect. 

The interpretation of the Wisconsin events’ 
which are also listed as supporting evidence for 
the particle is not convincing. In each of these 
events a K~ meson of 40 Mev is produced alone 
by a singly charged particle which may be at 
rest. The K mesons are arrested in the stack 
and they give rise to p endings. Yamanouchi 
proposes that they are due to an interaction of 
D~ mesons according to the scheme, 

7 0.200 

D om, a” * +K~. 
Such an interpretation does not explain the 
equality of the K-meson energy in each event 
because the momentum of the bound neutron 
would considerably broaden the spectrum; evap- 
oration prongs would also be expected. A better 
interpretation is that they are due to the decay 
of a heavy particle as originally suggested by the 
authors.’ 

There have been reports of other mass 1450 
particles; in 1955 Menon reported the result of 
Fowler and Perkins at the Rochester Conference.® 
Mass determinations by grain density and scat- 
tering measurements on flat particles from 
cosmic-ray stars gave some indication of a peak 
at 1480m,~. When these measurements became 
more accurate, however, in larger stacks where 
the particle ranges could be determined, these 
mass values were not confirmed. 

An event which can be interpreted as a D*- 
meson decay according to the scheme suggested 
by Yamanouchi has been observed by Sen Gupta 
and Sinha® in a cloud chamber. 

In conclusion we may say that although a num- 
ber of hitherto unexplained events can be in- 
terpreted if such a particle is postulated, con- 
crete evidence for its existence is lacking and 
is certainly not provided by the anomalous Kt 
decays. 

We are grateful to Professor H. K. Ticho, 

Dr. J. H. Mulvey, Professor J. Orear, and 
Professor E. C. G. Sudarshan for many dis- 
cussions. We also wish to thank Professor J. 
Orear for prepublication information on the 
mass value of the Columbia anomalous K* decay. 
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PARTIAL-WAVE DISPERSION RELATIONS FOR MESON-NUCLEON SCATTERING* 


Reinhard Oehme 
Enrico Fermi Institute for Nuclear Studies and Department of Physics, University of Chicago, Chicago, Illinois 
(Received February 3, 1960) 


It is of interest to compute the implications for 
the pion-nucleon scattering amplitude of possible 
resonance effects in the two-pion system.’ In 
this note we propose a form of the az -N partial- 
wave dispersion relations* which is suitable for 
this purpose. Our formulas can be easily gen- 
eralized for amplitudes involving other particles 
with unequal mass. 

We write the covariant 7 -N amplitude in the 
form’ F =A -iy-3(k+k’)B, where Ag, =A*bq 
+A~3[Tq, Tg), etc. and k+p=k’+p’, k®? =k” =-?, 
p? =p” =-m?. In the following we discuss ex- 
plicitly only the invariant amplitudes A*; the 
extension to B+ and to center-of-mass quantities 
is straightforward. Let us first consider the 
A+ as functions of the invariant variables 

=-(k+p)’, ¢=-(k-k’)?. We assume that they 
are analytic in both variables except for the 
absorptive singularities* due to the possible in- 
termediate states of the three physical reactions 
associated with the 7 -N Green’s function. Then 
we have the branch lines*® 


z=s2(m+pu)?, Z=52(m+y)’, 


and ¢=t2(2u)?, (1) 
where Z=2m*+2u?-z-¢. The single nucleon 
poles z =m’, Z=m? appear only in the amplitudes 
B*, The variables z and { may be expressed in 
terms of g (¢=c.m. momentum) and 6 (c.m. 
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angle) by 
z=m+p?+2@+2[(P+m*)(g+u7)!”, (& 
= -2q¢°(1 -cosé). (2b 


Note that the complex z plane is mapped intoa 
Riemann surface with two leaves which are 
connected through the cut -m?<q?<-u?. We 
define sheet I (II) by the requirement that the 
root in Eq. (2a) is positive (negative) for real 
q>-wu*. Inthe z plane, the exterior of the 
circle |z| =m? - u? corresponds to sheet I and 
the interior to sheet II. 

Let us now consider the amplitude A as an 
analytic function of ¢ for -1<cos@<+1. It is 
regular on the Riemann surface mentioned above 
except for cuts along the real axes in both sheets. 
In particular, the partial-wave amplitudes, 


+1 
A (@)=2 1 90089 P (c086)A(¢, cos), (3) 


have the following cuts: (7N) from 0 to +o in], 
due to s2(m+u)*; (7N) from 0 to +o in I and 
from -« in II up to -m? and then down to -« in 

I, both due to 52(m+u)*; (x7) from -p? to -« 
in Iand I, due to ¢2(2u)*. The weight functions 
associated with the unphysical branch lines (#N) 
and (77) can be expressed in terms of the absorp- 
tive parts N(a,8) and M(a,8) (a,8 being the 
squares of energy momentum transfer), cor- 
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responding to the reactions 7+N~7+N and 
1+8-N +N, respectively. In the physical sheet 
we can write partial-wave dispersion relations 
inthe simple form 


eo ImA (q’) 
1 Pee ind 
Af?) -i/ dq q” tf 7 
72 
1 -2 ‘ ImA (9 ) 
T bi dq q” -¢ , (4) 


which should be amended by the usual substrac - 
tions. However, in the interval -m?<q@< -y? 

the real weight ImA 7 can be expressed in terms 
of the absorptive parts N and M only indirectly 
by writing a similar dispersion formula in sheet 
I. Dispersion relations involving both sheets 
can be easily written down, but it seems to be 
advantageous to use the invariant variable z for 
those portions of the cuts which map onto the 
real axis in the z plane, and to retain the vari- 
able g for that part of the branch line (77), 

which maps onto the circle |z| =m? -?. Omit- 
ting again appropriate subtractions, we can 





c(s) 


2m? +2u?-s-5 





write 
oc ImA (s) 
= | J 
Afed=2f ads $-2 
1 (m-u)* ImA ,(s) 
rn s ere a na 
+F{e)+F, (e’), (5) 
where 
iH 2s(d) M {(s(a)) 
F (2) wale fa s(A) -m? -p?-2d s(A)-z ’ 
(6) 
with 
(A) = m? + up? +22 +2i[(A +m?)(-a - uw?) P, (7) 


and, for s<(m- uy)’, 


ImA 7 8) = 6(-s)e(s +m? - u2)N 7 (8) +M 7 (8): (8) 


The functions M7 and N.7 are given in terms of 
the absorptive parts M(a,8) and N(a, 8) by 


M {(s)=- 75 Sus dt P (1 +t/2q)M(t, s), 





N{s)=-gal, asp, lis 


m’+2u-s 


where ¢ = g*(s) is given by the inverse of Eq. (2a) 
with z=s; and c(s)=(m? - u?)?/s for 0,s<(m - zp)’, 
c(s)=(m+p)? for s$O_. Theterm Fy+Fy in 

Eq. (5) contains the low-energy properties of the 
4 system. Using Eqs. (6) and (9), it can be 

easily written in terms of real and imaginary 
parts of M(t, s).> The latter quantity is complex 
due to the appearance of s(A) given by Eq. (7). 

ha partial-wave expansion of M(t, s) the vari- 
able s will appear only in the argument of the 
spherical harmonics. 
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S. W. MacDowell, Phys. Rev. 116, 774 (1959). 

3See, for example, G. F. Chew, Annual Reviews of 
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1959), Vol. 9, p. 29; this paper contains further 
references. 

‘R. Oehme, Nuovo cimento 13, 778 (1959). 

‘The absorptive parts M(a,8) and N(a,8) can also 
become complex for large negative values of 8. The 
contributions from the corresponding imaginary parts 
cancel each other in Eq. (5). 
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CLASSICAL ELECTRODYNAMIC EQUATIONS OF MOTION WITH RADIATIVE REACTION* po 
Gilbert N. Plass of | 

Aeronutronic, Ford Motor Company, Newport Beach, California dia 

(Received November 5, 1959) the 


The classical relativistic equations of motion 
of a charged particle including the force of ra- 
diative reaction have been derived by Dirac’ and 
by Wheeler and Feynman.” The existence of 
divergent solutions of these equations which in- 
crease exponentially with time has been dis- 
cussed extensively in the literature.'»* The ex- 
istence of these solutions has been the basis for 
the fairly widespread belief that these equations 
provide only an approximate description of the 
behavior of a classical radiating charged par- 
ticle. Eliezer*® has claimed that no physical 
solutions for these equations exist for (1) the 
field of a thin infinite charged plate, (2) an 
attractive Coulomb force, and (3) a repulsive 
Coulomb force. On the contrary it can be shown 
that a physically reasonable solution for these 
equations exists for any force field that satisfies 
certain conditions given below. 

For simplicity consider first the one-dimen- 
sional, nonrelativistic equation of motion, 


ma-mb-‘a=f(t), (1) 
where 
b-+ =2e?/3mc5, (2) 


f(t) is the force acting on the particle, a is the 
acceleration of the particle, and dots signify 
differentiation with respect to time. 

The exact nondivergent solution of Eq. (1) is 


ma(t) =b of eo Ob 4s tat". (3) 
fv] 


Since Eq. (1) contains a third derivative with 
respect to the position of the particle, the appro- 
priate physical boundary conditions require that 
the initial acceleration as well as the initial 
velocity and position be specified. The initial 
value of the acceleration is obtained by setting 
t=0 in Eq. (3). It is seen from Eq. (3) that the 
acceleration at a given time is determined by the 
force that will act on the particle over a time 
interval of the order of b~* (which is approxi- 
mately 10-** sec) in the future. This phenomenon 
of pre-acceleration has been shown’’? not to 
violate our ordinary causal concepts over time 
intervals greater than b™’. 

From the theory of Laplace transforms it fol- 


248 


lows that a nondivergent, physical solution of yj} eq 
(1) for the position, velocity, and accelerationg) nw 
the particle exists whenever (1) f(#) is continuoy} ca: 
or piecewise continuous; (2) t’”|f(t+t’)| is (3) 
bounded near ¢’ =0 for some number ”, where ot! 
n<1; (3) e~ 5t’ |f(t+ t’)| is bounded for large va-| mi 
ues of t’. The solution for the three-dimension} ha 
nonrelativistic equation has the same form as pl 
Eq. (3) with a and f replaced by the en 
vector quantities. 

A general solution cannot be given for a space. 
dependent force. Whenever a solution exists of 
Eq. (3) for f[r(é)], it is nondivergent provided } us 
that the above conditions on the force are satis-} cl 
field. The existence of such a solution cannot} {fo 
proved in general; however, nondivergent solu- | T! 
tions have been obtained for a number of partic- 
ular space-dependent forces, some of which are 
listed below. 

An alternative solution to Eq. (1) is 


se B 


_ = ,-n,(n) 
ma(t)= 246 fr", (4 


where fw is the nth derivative of f. Thus, the 
acceleration at the time ¢ can also be obtained 
from all of the derivatives of the force evaluated 
at the same time ¢. This solution is usually very 
rapidly convergent, since a physical force sel- 
dom changes appreciably over a time interval 
iia 

Physically acceptable solutions exist for re- 
lativistic particles whenever a similar solution 
exists for the same force field at nonrelativistic 
velocities. The one-dimensional relativistic 
equation of motion is 


» te 3 oe» u?\¥2 
iF ity amt a) 


where u is the derivative of position with respect 
to the proper time, 7, and dots in this equation 
indicate differentiation with respect to 7. The 
exact, nondivergent solution of this equation is 


T 
u=Cc sinh{sinh (ug/c)+ (me) 1 f(r')dt’ 
0 





(me) tf eT rer) -f(r)]dr'}, 
0 
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where “, is the value of u at r=0. The first two 
terms on the right side of Eq. (6) are the solution 
of the relativistic equation of motion without ra- 
diative reaction, while the last term represents 
the effect on the motion of radiative reaction. 

The exact, nondivergent solutions of these 
equations of motion have been obtained for a 
number of particular force fields including the 
cases where (1) f=C(t-t,)”, n2-1; (2) f=C sinwt; 
(3) f=C exp|-(t-t,)?/207]; (4) f=C, x .<x<x,; f=0 
otherwise; (5) crossed, oscillating electric and 
magnetic fields; (5) constant magnetic field; (6) 
harmonic oscillator; (7) field of a thin charged 
plate. Furthermore, by numerical procedures, 
nondivergent solutions have been obtained for 
three-dimensional motion under the influence of 
attractive and repulsive Coulomb forces. 

Thus, when the proper boundary condition is 
used to specify the initial acceleration of a 
charged particle, nondivergent solutions exist 
for most force fields of interest in physics. 
Therefore, there no longer appears to be any 


reason for not accepting these equations of mo- 
tion. To the best of our knowledge, they provide 
an exact description of the motion of a charged 
particle within the framework of classical theory. 

The author wishes to thank Sylvan Katz for 
several valuable discussions during which the 
present method for obtaining the mathematical 
solutions for the equations of motion was for- 
mulated. Professor John Wheeler made many 
useful suggestions during the early stages of this 
work. 
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ABSTRACTS 








In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 


SIGNIFICANCE OF POTENTIALS IN QUANTUM 
THEORY. W. H. Furry and N. F. Ramsey, Ly- 
man Laboratory of Physics, Harvard University, 
Cambridge, Massachusetts (Received November 
16, 1959). 


The effects of the scalar and vector potentials 
in quantum mechanics, which were pointed out 
recently by Aharonov and Bohm, are discussed 
from the point of view of the consistency of the 
quantum- mechanical description of interference 
experiments. A well-known requirement for 
this consistency is that if any measuring device 
is introduced that can be used to determine 
which path the particle has taken, it must have 
the effect of eliminating the interference phenom- 
enon. Two conceptual experiments are discussed 
corresponding to the two phase effects noted by 
Aharonov and Bohm. In each case it is found 
that the phase effect is of just the magnitude re- 
quired to destroy the interference pattern when 
the circumstances are such that no pattern 
should be observed. 


RELATION BETWEEN INELASTIC NEUTRON 
SCATTERING AND THERMODYNAMIC FUNC- 
TIONS OF LIQUID HELIUM. Michael Cohen, 
University of Pennsylvania, Philadelphia, Penn- 
sylvania (Received October 26, 1959). 


A model of liquid helium is analyzed, in which 
the liquid is regarded as a collection of excita- 
tions (“rotons” only with energy >A) with an ar- 
bitrary pairwise number-conserving interaction. 
The entropy and normal fluid density of the li- 
quid and the energy distribution of scattered neu- 
trons are computed as power series in the den- 
sity of excitations exp(-A/kT). The first terms 
containing effects of the interactions are studied. 
When the interactions are weak, the entropy 
[through order exp(-24/kT)] is simply related to 
the neutron scattering, the connection being 
correctly given by the formula of Bendt, Cowan, 
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and Yarnell. For strong interactions there ap- 
pears to be no simple connection. Even when 
interactions are weak, the first correction to th 
normal fluid density involves information which 
is not contained in the neutron scattering. 

A method due to Bloch and de Dominicis is 
used in the analysis, and leads to a new form 
for the second virial coefficient. This is closely 
related to a curious new form for the level shift 
of a particle in a large spherical box, under the 
influence of a central potential. 


ELECTRIC FIELD DISTRIBUTIONS IN AN 
IONIZED GAS. I. Bernard Mozer* and Michel 
Baranger, Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania (Received November 
23, 1959). 


A method previously described is used to cal- 
culate the probability distribution of the low- 
frequency component of the electric field at a 
neutral point, the distribution of the low-fre- 
quency component at an ion, and that of the high- 
frequency component at an electron. The results 


are compared with those obtained by other author, 


*Now at Brookhaven National Laboratory, Upton, 
New York. 


SELF-CONSISTENT FIELD AND MOTION OF 
ELECTRONS WHICH HAVE A RANGE IN CA- 
NONICAL ANGULAR MOMENTUM IN A UNI- 
FORM MAGNETIC FIELD. Lewi Tonks, Law- 
rence Radiation Laboratory, University of Cali- 
fornia, Livermore, California (Received Noven- 
ber 2, 1959). 


The self-consistent theory of relativistic elec- 
trons circulating in a uniform impressed mag- 
netic field, the “Astron problem, ” has been 
generalized to the extent that a range of canonical 
angular momentum among monoenergetic elec- 
trons has been treated. For simplicity, the den- 
sity distribution in phase space has been chosen 
to be uniform over a finite momentum range. 
Just as in the single- momentum case, field re- 
versal is found, but new field and spatial density 
configurations appear. The uniform distribution 
is found to be consistent with isotropic regions 
of constant spatial density and constant magnetic 
field. The thickness of transition layer between 
vacuum and such a region conforms, within 
limits, to an empirical relation previously 
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found. The limit to the number of electrons per 
wit axial length of layer still exists. The curves 
relating the ratio of internal to external field to 
the layer strength still show multiple values of 
poth ratio and strength in certain ranges. Tra- 
jectories have been calculated and plotted for 
several cases. 


GROUND-STATE ENERGY OF A MANY-FERMION 
SYSTEM. W. Kohn, Carnegie Institute of Tech- 
nology, Pittsburgh, Pennsylvania, and J. M. 
Luttinger, University of Pennsylvania, Phila- 
delphia, Pennsylvania (Received October 19, 

1959). 


This paper contains a critique of the Brueckner - 
Goldstone perturbation series for the ground - 
state energy of an interacting gas of Fermions. 
We have calculated this energy by first construc- 
ting the grand partition function at finite tem- 
perature, and then carefully taking the limit as 
T-0. In general this leads to a series which 
differs from that of Brueckner and Goldstone. 

An exception is the case where both the unper- 
turbed single-particle energy as well as the 
interaction potential have spherical symmetry. 
Reasons for the breakdown of the Brueckner- 
Goldstone formalism are briefly discussed. 





QUANTUM STATISTICS OF NONIDEAL SYS- 
TEMS. Frank H. Stillinger, Jr.,* and John G. 
Kirkwood, Department of Chemistry, Yale 
University, New Haven, Connecticut (Received 
October 23, 1959). 


Anew cluster development for the logarithm 
of the grand partition function of a system of 
interacting particles is derived. The leading 
term in this expansion is the pressure exerted 
byan ideal Bose or Fermi gas at the same tem- 
perature and absolute activity Z as the actual 
system. Succeeding terms involve quantum 
(luster integrals which themselves depend upon 
Z, unlike their classical analogs. The definition 
ofthese cluster integrals follows in a natural 
fashion using techniques illustrated by construc- 
tion (in closed form) of the successive Z deriva- 
tives of the Bose and Fermi ideal gas grand 
fartition functions. It is not possible (except 
inthe classical limit) to eliminate Z explicitly 
tetween the pressure and density series, so 


that the equation of state must remain in para- 
metric form. 


*Present address: Bell Telephone Laboratories, 
Murray Hill, New Jersey. 
TDeceased. 


TRANSPORT PHENOMENA IN SLIGHTLY ION- 
IZED GASES: HIGH ELECTRIC FIELDS. 
Mahendra Singh Sodha, Physics Division, Ar- 
mour Research Foundation, Chicago, Illinois 
(Received November 9, 1959). 


Starting with the electron velocity distribution 
obtained by Chapman and Cowling for a Lorentz- 
ian gas, in the presence of an electric field, the 
author has investigated the variation with electric 
field of a number of transport properties, arising 
from a magnetic field, perpendicular to the 
electric field and temperature gradient in the 
gas. The applicability of the results to semi- 
conductors has also been pointed out. A constant 
mean free path has been assumed, which is 
validated by experiments for helium. 


LOW-LYING EXCITATIONS IN A BOSE GAS OF 
HARD SPHERES. F. Mohling, Department of 
Physics, Columbia University, New York, New 
York, and A. Sirlin, Department of Physics, 
New York University, New York, New York 
(Received November 2, 1959). 


The pseudopotential method is used to calcu- 
late the low-lying excitation energies of a Bose 
gas of hard spheres at T=0 to an order beyond 
that previously calculated by this technique. The 
results are in agreement with those obtained by 
Beliaev using a different approach. The phonon 
velocity is found to be equal to the velocity of 
compressional waves to the order of approxi- 
mation considered. 


STATISTICAL MECHANICAL THEORY OF 
FERROMAGNETISM: HIGH-DENSITY LIMIT. 
R. Brout, Laboratory for Atomic and Solid State 
Physics and Department of Physics, Cornell 
University, Ithaca, New York (Received Novem- 
ber 30, 1959). 


The partition function of the Ising model of 
ferromagnetism is examined in the limit of high 
density in the anticipation that in the limit of in- 
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finite density one recovers the Weiss molecular 
field. The formal parameter of expansion is 
1/z, where z is the number of spins in the range 
of the exchange potential (not restricted to near- 
est neighbor interactions). In the absence of 
long-range order, only ring diagrams in the 
cluster expansion contribute. These give a di- 
vergence in the specific heat at RT, = 2) sai Vii 
where vj; is the exchange potential. This is the 
molecular field value for the Curie point T;. In 
the presence of a magnetic field the partition 
function is evaluated for fixed magnetic moment 
M in the same approximation, M being deter- 
mined by minimization. This results in a sus- 
ceptibility differing from the molecular field 
theory and hence an inconsistency in the theory. 

The inconsistency is traced back to the obser- 
vation that the acceptance of ring diagrams is 
equivalent to the Gaussian model of Kac and 
Berlin which violates the sum rule Du?) =N. 
Hi; is the “spin” per particle and N is the total 
number of particles. This condition is reinstated 
by insuring the sum rule. The result leads to 
the spherical model. Thus, a consistent high- 
density approximation to the Ising model is the 
spherical model. Below the Curie point or for 
fixed magnetic field, M is again held fixed and 
only the Fourier components of the spin den- 
sity with nonvanishing wave vector are “spheri- 
calized.” The result leads to a physically ac- 
ceptable model which becomes the molecular 
field theory at low temperatures or high fields 
and deviates in O(1/z) in general. Formally, 
the results are simply expressed in terms of a 
temperature-dependent Weiss field. These re- 
sults differ from the ordinary spherical model 
which is physically unacceptable below the Curie 
point. However a molecular field modification 
of the spherical model due to Lax yields the 
same result when properly interpreted. 

It is shown that the above results are also 
valid (to the same approximation) in the quantum 
mechanical Heisenberg model, for temperatures 
above the Curie point. 


EXCITATION OF SPIN WAVES IN FERROMAG- 
NETS: EDDY CURRENT AND BOUNDARY 
CONDITION EFFECTS. P. Pincus, Department 
of Physics, University of California, Berkeley, 
California (Received November 16, 1959). 


The boundary condition for the transverse 
magnetization is derived when there is a surface 





252 


—, 


anisotropy field H,. Writing the transverse ma. 


netization as a sinkz +8 cos kz, we have 


— 
a/B =-tka+ (H ,/H ak), 


where H, is the exchange field and a the lattice 
constant. A similar result is found when there 
is an antiferromagnetic surface layer. For 
H,#0 spin-wave modes can be excited by a uni- 
form rf field in a ferromagnet. The power ab- 
sorbed in each mode in an insulator is calculate 
as a function of the surface anisotropy field. 
The excitation of the exchange modes is calcu- 
lated for a metal with eddy-current damping. 
The eddy currents are found to have a large 
effect only on long-wavelength spin waves. The 
line shape in a thick metal plate is calculated 
for H, normal to the plate. 


NUCLEAR QUADRUPOLE AND ELECTRONIC 
HEAT CAPACITIES OF BISMUTH. Norman E. 
Phillips, Department of Chemistry and Lawrence 
Radiation Laboratory, University of California, 
Berkeley, California (Received November 17, 
1959). 


The heat capacity of bismuth has been meas- 
ured from 0.1 to 2.0°K and found to be repre- 
sented by C =2.8x10-77~* +2.1 x10 °T +1944(7/6) 
joules /mole deg with 6=120.4-0.677. The T” 
term in the heat capacity is assumed to be asso- 
ciated with the alignment of the nuclear electric 
quadrupole moment in the electric field gradient 
of the crystal and is used to obtain the value 
25 Mc/sec for the quadrupole coupling constant. 
The linear term is used, together with known 
parameters for the electrons in the conduction 
band, to obtain an average effective mass for 
the holes in the valence band of 0.9 times the 
free electron mass. 


EFFECT OF OXIDATION ON THE CHARACTER- 
ISTIC LOSS SPECTRA OF ALUMINUM AND 
MAGNESIUM. C. J. Powell* and J. B. Swan,! 
Department of Physics, University of Western 
Australia, Nedlands, Western Australia (Re- 
ceived October 26, 1959). 


Measurements of the characteristic electron 
energy loss spectra of aluminum and magnesium 
were made (in a reflection experiment) during 
oxidation of a fresh evaporated layer of either 
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metal. It was found that surface oxidation re- 
sults in the rapid disappearance of the low-lying 
energy losses (10.3 ev in aluminum and 7.1 ev 

in magnesium) and the appearance of modified 
low-lying losses of 7.1 ev in aluminum and 4.9 
evin magnesium. The general changes in the 
loss spectra and the particular changes in the 
spectrum of aluminum were in good agreement 
with the predictions of Ferrell and Stern. 


*Now at the Imperial College of Science and Tech- 
nology, London, England. 

t Now on leave at Brookhaven National Laboratory, 
Upton, New York. 


LIMITING PROPERTIES OF NUMBERS OF SELF- 
AVOIDING WALKS. J. M. Hammersley, Institute 

of Statistics, Oxford, England (Received Novem- 

ber 16, 1959). 


A counter-example is given to a conjecture of 
Fisher and Sykes [Phys. Rev. 114, 45 (1959)] on 
the number of self-avoiding walks arising in the 
excluded-volume problem in the Ising model of 
ferromagnetism. 


NUCLEAR QUADRUPOLE SPIN-LATTICE RE- 
LAXATION IN ALKALI HALIDES. E. G. Wikner,* 
W. E. Blumberg,! and E. L. Hahn, Department 

of Physics, University of California, Berkeley, 
California (Received November 25, 1959). 


Nuclear quadrupole spin-lattice relaxation 
times have been measured in alkali halide crys- 
tals by the pulsed magnetic resonance technique. 
Measurements were made on Na”* in NaCl, NaBr, 
and NaI; Cl°° in NaCl and KCl; Br’®»® in NaBr, 
KBr, RbBr, and CsBr; Rb® in RbCl and RbBr; 
and I’? in NaI, KI, and CsI. Over a tempera- 
ture range of 298°K to 195°K the relaxation times 
are inversely proportional to the square of the 
absolute temperature. The data are compared 
to relaxation times calculated from an ionic 
crystal model of Van Kranendonk and a covalent 
model of Yosida and Moriya. The ionic model is 
modified to include the interaction between the 
tuclear quadrupole moment and the electric field 
gradient due to electric dipole moments asso- 
ciated with optical modes of vibration. Neither 
of these models alone predicts the experimental 
telaxation times for all cases, but a combination 
of the two effects is required. The modified ionic 





model applies reasonably well to crystals which 
contain the lighter ions. 

‘Present address: General Atomics, San Diego, 
California. 


TPresent address: Bell Telephone Laboratories, 
Murray Hill, New Jersey. 


CRYSTAL STRUCTURE OF FERROELECTRIC 
LiH,(SeO,),. K. Vedam, Y. Okaya, and R. Pepin- 
sky, Department of Physics, Crystal Research 
Laboratory, The Pennsylvania State University, 
University Park, Pennsylvania (Received Novem- 
ber 30, 1959). 


The structure of the room-temperature ferro- 
electric LiH,(SeO,), has been determined by 
X rays, using the heavy-atom method, and re- 
fined on the IBM 704. The crystals are mono- 
clinic, with space group Pn and a=6.25, A, 
b=7.88, A, c=5.43, A, 8=105.2°. Fairly strong 
O-H:--O bonds with distances 2.52, 2.56, and 
2.57 A are found, nearly perpendicular to the 
polar direction. The O-Se-O angles in one of 
the two selenite ions are rather similar; in the 
other ion these angles are unequal, as in the 
structure of H,SeO,. Possible positions for the 
Li ions are given based on crystal-chemical 
considerations. 


SPIN RESONANCE OF CHARGE CARRIERS IN 
GRAPHITE. G. Wagoner, Research Laboratory, 
National Carbon Company, Cleveland, Ohio (Re- 
ceived November 25, 1959). 


The observations reported here of the electron 
spin resonance in quite perfect single crystals 
of graphite clearly establish that the resonance 
arises from mobile charge carriers. The line 
shape is of the Dyson form which is character - 
istic of conduction electron spin resonance in 
metals. The intensity of the spin resonance 
agrees, both in absolute magnitude and in tem- 
perature dependence, with values calculated 
from the band model of graphite by McClure. 
The g value of the resonance shows a remark- 
ably large anisotropy which depends strongly on 
temperature and on the position of the Fermi 
level with respect to the band edge. At room 
temperature in pure graphite, g varies from 
2.0026 + 0.0002 to 2.0495+0.0002 as the magnetic 
field is shifted from perpendicular to parallel to 
the c axis. The g-value anisotropy increases 
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with decreasing temperature; g, becomes 2.127 
at 77°K while g, remains constant. The line 
width of the resonance is a few gauss (T, = 2.0 
x107® sec) which is extremely narrow in com- 
parison with the field shifts caused by changes 
of anisotropy with temperature. This indicates 
that for conduction states in graphite, the g 
value is a strong function of the wave vector and 
that the line is narrowed by an averaging pro- 
cess ink space. This averaging is similar to 
that which occurs in motional and exchange 
narrowing. 


NORMAL MODES OF A LATTICE OF OSCILLA- 
TORS WITH MANY RESONANCES AND DIPOLAR 
COUPLING. U. Fano, National Bureau of Stand- 
ards, Washington, D. C. (Received November 

2, 1959). 


The normal modes of a lattice of coupled di- 
poles are studied as a model of the collective 
excitations of electrons in condensed materials. 
Two types of oscillations are found, in which 
electrostatic coupling has a dominant influence. 
One of them is analogous to the oscillation of an 
electron plasma and has a high dipole moment. 
Other collective oscillations have a low net di- 
pole moment, owing to destructive interference 
between out-of-phase components. These two 
types of oscillation occur in systems with a 
sufficiently high density of oscillator strength 
in space and in spectrum. A simple estimate 
indicates that most condensed materials fulfill 
this condition. 


PULSED NUCLEAR RESONANCE SPECTROS- 
COPY. M. Emshwiller,* E. L. Hahn, and D. 
Kaplan,! Department of Physics, University of 
California, Berkeley, California (Received 
November 2, 1959). 


A technique for the detection of weak nuclear 
resonance interactions in solids is carried out 
by a pulsed method which obtains nuclear double 
resonance. The resonance of the unknown spe- 
cies to be detected does not require a Boltzmann 
population difference in spin orientation, but 
must have a sufficient dipole-dipole interaction 
with a second spin species. At resonance, a 
single 180° pulse reorients the unknown spins at 
the time of the 180° pulse in the 90°-180° pulse 
sequence necessary to obtain the observed spin- 
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echo signal of the second species. A reduction J jon 
of the spin-echo signal signifies double reso- ins 
nance due to changes in local dipolar fields, the 
coupled to the observed spins, which scrambles [ the 
their precessional phases. Nuclear quadrupole f ma 
coupling interactions of K, Cs, and Rb isotopes pho 


are measured in the chlorates of these ions, ger 
where the Cl** nucleus provides the observed cri 
nuclear quadrupole echo. An analysis is pre- ene 


sented for the case of low concentration of un- pri 
known spin species. Double quantum transitions } fr 
and special properties of nuclear quadrupole sili 


spectra are observed. pho 
* is j 
Present address: Bell Telephone Laboratories, 

: dep 
Whippany, New Jersey. ha 
tPresent address: Lockheed Missile Systems Divi- . 

sion, Palo Alto, California. 
TH 


INTERNAL FIELD EMISSION AT NARROW PE 
SILICON AND GERMANIUM p - n JUNCTIONS. Phy 
A. G. Chynoweth, W. L. Feldmann, C. A. Lee, § £a 
R. A. Logan, and G. L. Pearson, Bell Telephone f 1% 


Laboratories, Murray Hill, New Jersey, and T 
P. Aigrain, Ecole Normale Supérieure, Paris, per 
France (Received November 5, 1959). cd 


A detailed study has been made of the reverse § 420 
characteristics of several silicon and germaniunf ‘i 
alloyed p - n junctions with breakdown voltages ff “ 
in the range of about 0.1 to 0.8 volt. In these COp 
junctions the reverse current is generated al- §— ™ 
most entirely by internal field emission (tunnel- — ™' 
ing). The reverse bias characteristics are in- T 
sensitive to the dislocation density present so del 
that the tunneling current occurs mainly in un- f % 
distorted material. From capacitance studies tal 
it is established that these narrow junctions are ma 
very close to being ideal step junctions. The 


room temperature reverse characteristics are ten 
analyzed in terms of the usual tunneling pro- ne 





bability expressions and in particular, good 
agreement, both qualitative and quantitative, is 
found between experiment and theory. The tun- 
neling probability exp(-ae*?/E), when compared 
with experiment, yields values for ae® in agree 
ment with the theoretical ones to within a factor 
of less than 2 for both silicon and germanium. 
The critical voltage (the reverse bias voltage 
necessary to maintain a constant tunneling cur- 
rent) was measured as a function of temperature 
from 4.2°K up to temperatures as high as 700%. 
In germanium, the critical voltage drops mono 
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tnically as the temperature increases whereas 
insilicon, there is considerable structure in 
the curve. This is shown to be consistent with 
the tunneling being by direct transitions in ger- 
manium and by indirect transitions (involving 
phonon emission and absorption) in silicon. In 
germanium, the temperature dependence of the 
critical voltage arises from that of the direct 
energy gap while in silicon, it is determined, 
primarily, by the available phonon density. 
From an analysis of the temperature data for 
silicon that invokes the transverse acoustic 
shonons, the estimate of ae*? that is obtained 
is in excellent agreement with that found, in- 
dependently, from the analysis of the reverse 
characteristics. 


THERMOELECTRIC POWER OF DILUTE COP- 
PER ALLOYS. F. J. Blatt and R. H. Kropschot,* 
Physics Department, Michigan State University, 
fast Lansing, Michigan (Received October 12, 
1959). 


The thermoelectric power of dilute binary cop- 
per alloys containing one atomic percent Zn, Ge, 
Cd, In, Sn, and Sb was determined over the 
temperature range extending from about 8°K to 
320°K. The thermoelectric powers were ob- 
tained by measuring the thermoelectric emf’s of 
thermocouples formed of the alloys and pure 
copper. The absolute thermoelectric power of 
pure copper was obtained from measurements 
ona pure copper versus lead thermocouple. 

The results are analyzed in terms of the Frie- 
del theory. It is found that above 40°K satis- 
factory agreement with that theory can be ob- 
tained only if it is assumed that phonon drag 
makes a significant contribution to the thermo- 
electric power in the pure material below room 
temperature. The magnitude and temperature 
dependence of this assumed phonon drag contri- 
bution are in satisfactory agreement with theo- 
retical estimates. Moreover, the effect of alloy- 
ing on this phonon drag contribution also agrees 
vith theoretical estimates based on reasonable 
models. In the temperature range below about 
{0°K we have not been able to give a satisfactory 
interpretation of our results. Our measurements 
show that near 40°K the absolute thermoelectric 
bower of pure copper reverses sign, becoming 
legative, and attains an anomalously low mini- 
tum of about -1.8 microvolts/degree near 10°K. 
I'this behavior is due to the presence of minute 


amounts of impurities in the pure copper which 
give rise to a thermoelectric anomaly associated 
with the appearance of a resistivity minimum, 
then our results can be interpreted without much 
difficulty. However, the residual resistivity of 
the pure copper was so low that we do not believe 
that the low-temperature thermoelectric anomaly 
can be attributed to the presence of impurities. 

A research program on thermoelectric proper- 
ties of dilute alloys which should shed further 
light on these questions is outlined. 


"Present Address: National Bureau of Standards, 
Cryogenics Laboratory, Boulder, Colorado. 


LOW-TEMPERATURE IMPURITY CONDUCTION 
IN n-TYPE SILICON. K. R. Atkins, R. Donovan,* 
and R. H. Walmsley, Department of Physics, 
University of Pennsylvania, Philadelphia, Penn- 
sylvania (Received November 12, 1959). 


Hall coefficients and electrical resistivities 
were measured down to liquid helium tempera- 
tures for silicon specimens containing about 10*’ 
phosphorus impurities per cm® and about 10** 
boron impurities per cm’. The density of minor- 
ity impurities was determined during the pre- 
paration of the ingots, rather than deduced from 
the electrical measurements themselves. The 
results are extremely sensitive to the density 
of minority impurities. They are discussed in 
relationship to the theories of Conwell, Mott, 
and Price. 

= 


Now at Molecular Engineering Section, Westinghouse 
Air Arm, Blatimore, Maryland. 


LOW-TEMPERATURE INTERNAL FRICTION IN 
FACE-CENTERED CUBIC AND BODY-CENTERED 
CUBIC METALS. L. J. Bruner, * Department of 
Physics and Institute for the Study of Metals, 
University of Chicago, Chicago, Illinois (Re- 
ceived November 2, 1959). 


Data on the anelasticity produced by plastic 
deformation in various face-centered cubic and 
body-centered cubic pure metals and alloys are 
reported. Face-centered cubic materials stud- 
ied at temperatures from 4.2°K to 300°K include 
Cu, Al, and Al - 0.25 at. @ Cu. Body-centered 
cubic systems are Fe, Nb, and g brass. Bordoni 
peaks are observed in Cu and Al in agreement 
with-previous work, but are not found in either 
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strain aged Al-Cu alloys or pure Fe. A peak 
observed in Nb at 173°K is not believed to be a 
Bordoni type. Unexplained low-temperature 
internal friction peaks are also observed in 8 
brass. A new mechanism is proposed for dis- 
location relaxation in which the essential feature 
is the thermally activated motion of paired 
partial dislocations between vacancy -pinning 
points. It is in general qualitative accord with 
experiment, and permits semiquantitative evalu- 
ation of all essential parameters. 

“Now at International Business Machines Watson 


Laboratory at Columbia University, New York, New 
York. 


PROPERTIES OF p-TYPE InSb IN PULSED 
HIGH ELECTRIC FIELDS. M. C. Steele and M. 
Glicksman, RCA Laboratories, Princeton, New 
Jersey (Received November 6, 1959). 


The results of high electric field experiments 
on p-type InSb at 77°K are described. It is 
shown that electron-hole pair creation occurs at 
electric fields greater than 700 volts/cm. When 
a sufficient number of pairs are created the Hall 
coefficient changes from positive to negative. 
The question of whether holes or possible in- 
jected electrons initiate the pair creation is ex- 
amined in detail. An incipient negative resistance 
effect in transverse magnetic fields and the ab- 


sence of any self-pinch effects are also discussed. 


DIFFUSION AND IONIC CONDUCTIVITY IN 
CESIUM BROMIDE AND CESIUM IODIDE. David 
W. Lynch,* Department of Physics, University of 
Illinois, Urbana, Illinois (Received November 
9, 1959). 


The diffusion coefficients of the constituent 
ions were measured in CsBr and CsI single crys- 
tals between 300°C and 550°C by means of radio- 
active tracers and were compared with electrical 
conductivity measurements. Approximate satis- 
faction of the Nernst-Einstein equation indicates 
that the conductivity is nearly completely ionic 
and the diffusion measurements show that the 


halogen ion defects are the more mobile. Attempts 


to identify the mechanisms for ionic transport 
by means of the Bardeen-Herring correlation 
factor are discussed. Schottky defects seem 
likely but the assumption of an additional mech- 
anism for cationic transport is required. If 
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Schottky defects are predominant, their formatig, 
enthalpies are 2.0 ev and 1.9 ev in CsBr and Csl, 
respectively, and the activation enthalpies for 
halogen vacancy motion are 0.27 ev and 0.3 ev, 
respectively, while the cesium vacancy activatioy 
enthalpy is about 0.58 ev for both salts. 


"Present address: Institute for Atomic Research 
and Department of Physics, Iowa State University, 
Ames, Iowa. 


TEMPERATURE DEPENDENCE OF EDGE 
EMISSION IN CADMIUM SULFIDE. D. C. Rey- 
nolds, Aeronautical Research Laboratory, 
Wright-Patterson Air Force Base, Ohio (Re- 
ceived November 12, 1959). 


The temperature dependence of the blue and 
the green emission in cadmium sulfide has been 
investigated in the temperature range from 77K 
to 300°K. The blue emission shifts to longer 
wavelengths as the temperature is increased 
above 77°K with a temperature coefficient of 
4.9x107* ev/°C. The intensity of the green 
emission falls off exponentially above 77°K with 
very little shift in the spectral position of the 
green peaks. 


SPIN-PHONON INTERACTION IN PARAMAG- 
NETIC CRYSTALS. R. D. Mattuck* and M. W. P. 
Strandberg, Department of Physics and Research 
Laboratory of Electronics, Massachusetts In- 
stitute of Technology, Cambridge, Massachusetts 
(Received November 16, 1959). 


A general theory of the spin-phonon interaction, 
which is applicable to any iron group spin not in 
an S state, is developed. The theory employs 
a perturbation treatment that has a more direct 
physical meaning than techniques previously 
used and that leads to more accurate results. 
These results are presented in the form of an 
equivalent spin-phonon interaction Hamiltonian 
involving sums over products of spin operators 
and phonon creation-annihilation operators. The 
interaction between any two spin levels can then 
be calculated by using the spin wave functions 
associated with the usual “spin Hamiltonian.” 

It is shown that, owing to the dominant role 
played by the quadratic term in the above in- 
teraction, odd half-integer iron group spins 

(S> 4) obey quadrupole selection rules. A for- 
mula is derived for order -of -magnitude cal- 
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culations of the interaction strength. It is shown 
that acoustic experiments should provide the 
ideal way to test this theory in detail, and two 
methods of checking the quadrupole rule are 
proposed. Experimental results are reported 

on observed acoustic saturation in MgO doped 
with Cr**+*, on the absence of saturation be- 
tween low-field Kramers doublets in ruby, and 
an apparent saturation effect in F-center quartz. 


*Now at Royal Institute of Technology, Stockholm, 
Sweden. 


GALVANOMAGNETIC EFFECTS IN n-TYPE 
SILICON. W. E. Krag, Lincoln Laboratory, 
Massachusetts Institute of Technology, Lexing- 
ton, Massachusetts (Received November 3, 
1959). 


The galvanomagnetic effects in n-type silicon 
have been investigated experimentally and theo- 
retically at 77°K and 300°K. The samples ranged 
in resistivity from 0.2 to 27 ohm-cm (at room 
temperature). Measurements of the field de- 
pendence from 0.5 to 22 kilogauss of the Hall 
coefficient and the magnetoresistance for sev- 
eral orientations are presented, as well as 
measurements of the angular dependence of the 
magnetoresistance and planar Hall effects at 
similar fields. A number of theoretical calcula - 
tions using a constant relaxation time in the 
manner of Gold and Roth are presented which 
show the effect on the coefficients of a change 
in the mass anisotropy K,,. The Boltzmann 
transport equation, and collision times of the 
form 


vy =1 6 + LEM +N, 


v =1,6 + LE +N, 


where v7 and v; are the longitudinal and trans - 
verse collision times, Jj, J;, L, and N are con- 
stants, and & is the energy, have been used to 
calculate the galvanomagnetic coefficients for 
different amounts of lattice, ionized-impurity, 
and neutral-impurity scattering, and for dif- 
ferent anisotropies of the ionized-impurity scat - 
tering. The various components of the collision 
frequency were taken from available theoretical 
work. Detailed comparisons of the experimental 
data and the theoretical calculations are pre- 
sented which show that the assumption of a con- 
stant collision frequency allows one to predict 
qualitatively all the field dependencies and sym- 


metries which are found experimentally. It is 
also shown that with an anisotropic, energy- 
dependent collision frequency the quantitative 
agreement with experiment is considerably im- 
proved, especially with regard to the behavior 
of samples of different resistivities. 


SPIN RELAXATION AND LINE WIDTH IN AL- 
KALI METAL VAPORS. Arnold L. Bloom, 
Varian Associates, Palo Alto, California (Re- 
ceived November 16, 1959). 


Relaxation constants T, and T, have been com- 
puted for experiments involving optical pumping 
and optical detection in alkali metal vapor. The 
calculations have been performed for several 
possible spin relaxation mechanisms; namely, 
magnetic dipole, electric quadrupole, spin ex- 
change, and the optical pumping process itself. 
For all of these mechanisms a reorientation ex- 
periment will approximately predict a spin re- 
sonance linewidth (equivalent to the statement 
T,=T, for spin-1/2 particles). However, a spin 
reorientation experiment of the type originally 
performed by Dehmelt, employing circularly 
polarized light, gives ambiguous results because 
of the nonexponential character of the reorienta- 
tion. A more suitable experiment is one em- 
ploying hyperfine population differences and un- 
polarized light. 


CALCULATION OF THE MAGNETIC HYPER- 
FINE STRUCTURE COUPLING CONSTANTS OF 
NO. Héléne Lefebvre-Brion and C. M. Moser, 
Centre de Mécanique Ondulatoire Appliquée, 
Paris, France (Received November 18, 1959; 
revised manuscript received January 21, 1960). 


A calculation of the magnetic hyperfine struc - 
ture (hfs) coupling constants of NO has been 
made using an ab initio configuration interaction 
wave function built from LCAO-MO-SCF orbitals. 
The constants predicted by this wave function 
are in good agreement with experiment for the 
two constants related to the quantities (L2/7*) ay 
(constant a) and (sin?6S~/r*),, (constant d). 

The direct calculation of ((3 cos?6 -1)S,/7*) ay 
(constant c) gives a value which is about 25% 
less than that predicted by an approximation 
which implies that unpaired / and s are asso- 
ciated with the same electron. From this cal- 
culated value of c the quantity (y7(0)),, -3¢ 
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(constant 5) now is found to be about 40 Mc/sec 
instead of 68.91 Mc/sec which has been pre- 
viously used by experimentalists. 

The same wave function has been used to cal- 
cuiate the fine structure constant (spin-orbit 
coupling) for the 71 ground state of NO. The 
agreement with experiment is quite satisfactory. 

The calculation of the nuclear quadrupole coup- 
ling constant g gives about 73% of the observed 
electronic part assuming that @(N™) = 0.02 x10 


cm’. 


ELECTRON-HYDROGEN SCATTERING AT LOW 
ENERGIES. Takashi Ohmura* and Haruko 
Ohmura, * Division of Pure Physics, National 
Research Council, Ottawa, Canada (Received 
October 23, 1959). 


The effective range in the singlet system has 
been evaluated as 2.646+ 0.004 atomic units by 
using the asymptotic amplitude of the 202- 
parameter H wave function of Pekeris. This 
value of the effective range, together with the 
value of the electron affinity of H”, determines 
the scattering length in the singlet system as 
6.167. The effective range in the triplet system 
is calculated to be 1.219 atomic units by a 
Hartree-Fock approximation. It is shown that 
the effective-range approximation is very good 
for all energies at which only elastic scattering 
is allowed. The photo-ionization of H™ is briefly 
discussed on the basis of the effective-range 
theory. 


*on leave of absence from Department of Physics, 
University of Tokyo, Tokyo, Japan. 


RELATIVISTIC SELF-CONSISTENT SOLUTIONS 
FOR ATOMS OF LARGE ATOMIC NUMBER. 
Stanley Cohen, The RAND Corporation, Santa 
Monica, California, and Lawrence Radiation 
Laboratory, University of California, Berkeley, 
California (Received November 2, 1959). 


Relativistic self-consistent solutions, without 
exchange, have been obtained for several atoms 
of large atomic number by use of a general pro- 
gram for a high-speed computing machine. A 
short description of this program and of the self- 
consistent calculation is given. Eigenvalues for 
the individual electron subshells of the self- 
consistent mercury, tungsten, platinum, and 
uranium atoms are presented. A comparison of 
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the calculation with previous results for the 
mercury atom is also included. 


MAGNETIC HYPERFINE STRUCTURE OF THE 

GROUND STATE OF LITHIUM. R. K. Nesbet,* 

Lincoln Laboratory, Massachusetts Institute of 

Technology, Lexington, Massachusetts (Receive; 
November 30, 1959). 


The magnetic hyperfine splitting of the 7S 
ground state of the lithium atom is calculated. It 
is shown that the discrepancy between experi- 
ment and the value calculated in the traditional 
Hartree-Fock approximation can be accounted 
for quantitatively by the exchange polarization 
effect, which distorts one 1s orbital relative to 
the other. The present calculation obtains a 
value within one percent of the experimental 
value. A general procedure is proposed for 
evaluating operators that do not commute with 
the Hamiltonian, when approximate variational 
methods must be used. 


* 
Present address: Department of Physics, Boston 
University, Boston, Massachusetts. 


POSITRON ANNIHILATION IN AQUEOUS SOLU- 
TIONS. Georg Trumpy,* Institutt for Atomenergi, 
Lillestrgm, Norway (Received November 30, 
1959). 


The angular correlation of 2-quantum emission 
from the annihilation of positrons in different 
materials has been measured in an apparatus 
with 8 photon counters providing coincidences 
for 16 output channels. As positron targets were 
chosen indium, water, and aqueous solutions of 5 
paramagnetic salts and 10 other substances. It 
was confirmed that the amount of singlet posi- 
tronium formed is influenced by two processes: 
a reduction of positronium due to electron cap- 
ture by oxidizing substances and an increase of 
the triplet —- singlet conversion due to the elec- 
tron exchange with paramagnetic ions. 

The oxidation potential of positronium is found 
to be very nearly zero. The conversion rate 
seems to be proportional to the number of un- 
paired electrons on the dissolved ions. A dis- 
crepancy with the interpretation of Green and 
Bell for their lifetime experiments is discussed. 


*Now at A. B, Atomenergi, Studsvik, Tystberga, 
Sweden. 
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s{ELL-MODEL ASSIGNMENTS FOR THE 
fNERGY LEVELS OF Ct AND N*. E. K. War- 
purton and W. T. Pinkston,* Palmer Physical 
laboratory, Princeton University, Princeton, 
New Jersey (Received November 24, 1959). 


Electromagnetic transition widths, -reduced 
yidths, and inelastic scattering cross sections 
are calculated for the following states of N“ 
and C**; (1) the levels arising from the ground- 
state configuration, s*p*®; (2) the odd-parity 
levels arising from excitation of a 1p nucleon 
into the degenerate 2s,,. and 1d,,. shells; (3) the 
even-parity group of levels formed by excitation 
of two 1p nucleons into the 2s and 1d shells. The 
calculations for the s*p’° configuration are car- 
ried out using the wave functions of Elliott and 
of Visscher and Ferrell, and in jj coupling. 

The calculations for the odd-parity levels are 
done in the jj-coupling scheme. For the even- 
garity excited configuration an inert C’ core 
isassumed and M1 radiative widths are cal- 
culated for states arising from s?+d*+sd. The 
calculations are compared to the existing data. 
On the basis of this comparison, shell-model 
assignments are proposed for 19 of the 27 known 
levels below 11-Mev excitation in N“* and for all 
the known levels in C™ below 9-Mev excitation. 


‘Present address: Department of Physics and 
Astronomy, Vanderbilt University, Nashville, Tennes- 
see. 


ENERGY LEVEL PARAMETERS FROM NU- 
CLEAR RESONANCE FLUORESCENCE AT 7 
Mev. K. Reibel and A. K. Mann, University of 
Pennsylvania, Philadelphia, Pennsylvania (Re- 
ceived November 17, 1959). 


The recoil-broadened photon spectrum from 
the reaction F?°( p, ay)O** has been used to 
measure the elastic photon scattering cross 
sections at 7 Mev of 31 elements. The observed 
angular distributions are consistent with dipole 
transitions. A plot of the cross sections versus 
mass number shows definite peaks around the 
closed shell regions near Z =50, N=82 (Sn, Te, 
ad Ba), and Z =82, N=126 (Pb and Bi). For 

six medium and heavy elements self-absorption 
measurements were made which, when analyzed 
interms of a number of nonoverlapping Breit - 
Wigner resonances, yield values of the average 
jartial radiation widths to the ground states, the 












werage total radiation widths, and the average 
level spacings for those elements. The radiation 








widths are significantly larger than those de- 
termined from slow-neutron scattering and 
capture experiments and, excepting Pb and Bi, 
the average level spacings are also appreciably 
greater than would be expected from the neutron 
data. The observed widths and spacings are in 
order -of-magnitude agreement with the recent 
interpretation of the modified single-particle 
calculation of Blatt and Weisskopf. 


GAMMA-RAY INTENSITIES IN THE THORIUM 
ACTIVE DEPOSIT. Guy T. Emery* and Walter R. 
Kane, * Lyman Laboratory of Physics, Harvard 
University, Cambridge, Massachusetts (Received 
November 24, 1959). 


The relative intensities of the gamma rays of 
the descendants of radiothorium have been studied 
by measuring the photoelectron spectra from 
thorium and platinum foils. The relative sensi- 
tivity of the arrangement at different energies 
was found by using other sources having gamma 
rays of known relative intensity. Comparison of 
the measured gamma-ray intensities with the 
measurements by others of internal conversion 
intensities allows internal conversion coefficients 
to be computed and the multipolarities of several 
transitions to be determined. An electric mono- 
pole transition was found in Po**. The gamma- 
ray intensities are used to find the intensities of 
beta-ray branches. Gamma-ray intensities are 
compared with the known intensities of long- 
range alpha particles from Po*!” and transition 
probabilities are estimated for some electromag- 
netic transitions between states of that nucleus. 
The level schemes of Po?” and Pb? are dis- 
cussed in the light of the information found here 
and of other recent information. 


*Now at Brookhaven National Laboratory, Upton, 
New York. 


SLOW-NEUTRON SCATTERING BY THE TITA- 
NIUM ISOTOPES. C. G. Shull* and M. K. Wilkin- 
son, Oak Ridge National Laboratory, Oak Ridge, 
Tennessee, and M. H. Mueller, Argonne Na- 
tional Laboratory, Lemont, Illinois (Received 
November 17, 1959). 


Neutron diffraction studies are reported on 
isotopically enriched samples of TiO, from which 
are evaluated the coherent scattering amplitudes 
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of the titanium isotopes. Scattering amplitudes 
of +0.48, +0.33, -0.58, +0.08, and +0.55x10"" 
cm were established for the titanium isotopes of 
mass 46, 47, 48, 49, and 50, respectively. The 
major isotope Ti*® is thus responsible for the 
anomalous scattering amplitude, -0.34x10"™ 
cm, characteristic of the normal element. Pro- 
nounced nuclear scattering resonance effects on 
the observed neutron scattering are suggested 
to occur for the most of the isotopes. 


"Present address: Massachusetts Institute of Tech- 
nology, Cambridge, Massachusetts. 


SLOW-NEUTRON RESONANCE SPECTRO- 
SCOPY. I. U***. J. L. Rosen, J. S. Desjardins, 
J. Rainwater, and W. W. Havens, Jr., Depart- 
ment of Physics, Columbia University, New 
York, New York (Received November 20, 1959). 


The results of time-of-flight measurements of 
U*** resonances in the region 90 - 1300 ev are 
presented and resonance parameters for levels 
up to 1000 ev are obtained. Neutron widths for 
the 55 observed levels and radiation widths for 
32 of the stronger levels are deduced. The re- 
duced neutron width distribution is found to be 
in good agreement with the theoretical prediction 
of Porter and Thomas for a single channel pro- 
cess, while the level spacing distribution agrees 
with the “repulsion” formula suggested by Wig- 
ner. The average value of the radiation widths 
was found to be (24.6+0.8) x10-° ev, while the 
average reduced neutron width and level spacing 
were found to be (1.76 + 0.26) x 107° ev and 18.5 
+1.3 ev, respectively. These values are in good 
agreement with earlier results reported by other 
workers. A strength function of (0.95+ 0.15) x10~* 
is obtained. 

It appears on the basis of their size and num- 
ber, that several of the weaker levels may be 
due to p-wave neutrons. 


PRECISION MEASUREMENT OF THE TOTAL 
NEUTRON CROSS SECTION OF U*** BETWEEN 
0.000818 AND 0.0818 ev. G. J. Safford, W. W. 
Havens, Jr., and B. M. Rustad, Columbia Uni- 
versity, New York, New York, and Brookhaven 
National Laboratory, Upton, New York (Received 
November 20, 1959). 


The absolute value of the total neutron cross 
section of U*** was measured at neutron energies 
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between 0.000818 ev and 0.0818 ev for two typg 
of samples, a metallic foil and D,O solutions 
of uranium nitrate. Balanced solutions of 
U**°0,(NO,), and U**O,(NO,), were used to de- 
termine the difference between the total cross 
sections of U*** and U***. This value, when cop 
bined with the relatively small known value of 
the total cross section for U**, gives o7(U™) 
=587+5 barns at 0.0253 ev. The measurement 
on the metallic U?** foil agreed with the meas- 
ured total cross section determined from the 
liquid solution data to better than 1%, yielding 
o(U**) = 586 +2 barns at 0.0253 ev. 


SLOW -NEUTRON TOTAL AND FISSION CROs 
SECTIONS OF U***. M. S. Moore, L. G. Miller, 
and O. D. Simpson, Atomic Energy Division, 
Phillips Petroleum Company, Idaho Falls, Ida 
(Received November 25, 1959). 


The slow-neutron total and fission cross sec- 
tions of U** have been measured from 0.02 ev 


OL 


is 


10\ 
TIS 
Ene 


to 200 ev on the MTR (Materials Testing Reactuff, 


fast chopper. The strong resonances are re- 
solved below a neutron energy of 15 ev, and 
show marked interference effects in the fission 
cross section. No resonances are observed in 
the total cross section which are not also pres- 
ent in the fission cross section, except for thos 
attributed to the known contaminants in the 
samples. An estimate of the neutron strength 
function (T,,°/D), made by an area analysis, 
gives the value (1.0+ 0.2) x10™ for this energy 
region in U**’, 


MULTILEVEL ANALYSIS OF THE SLOW- 
NEUTRON CROSS SECTIONS OF U***. M. S. 
Moore and C. W. Reich, Atomic Energy Divi- 
sion, Phillips Petroleum Company, Idaho Falls, 
Idaho (Received November 25, 1959). 


A multilevel analysis of the slow-neutron total 
and fission cross sections of U?** from thermal 


energies to 11 ev has been carried out, under tit 


assumption that in one of the two spin states the 
nucleus undergoes fission primarily through 2 
small number of channels. The cases of one asi 
of two fission channels for this spin state have 
been investigated in detail. Although the analy- 
sis does not yield an exact value for the number 
of available fission channels, the data are ade- 
quately fit with a two-fission-channel formula. 
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jg not necessary to invoke any negative-energy 
sonances in the analysis. The presence of a 


waving a 1/v energy variation, is indicated. At 
nermal energies, this component accounts for 


sections. An interpretation of this component 
being due to one spin state of the compound 

weleus is presented. Under this assumption, 

ne fission characteristics of the two spin states 

of the compound nucleus formed by s-wave neu- 

rons are quite different. 


LOW-ENERGY NEUTRON CROSS SECTIONS OF 
FSSIONABLE NUCLEI. Erich Vogt, Atomic 
nergy of Canada Limited, Chalk River, Ontario, 
(anda (Received November 30, 1959). 


The method of analysis developed in a previous 
maper is applied to the low-energy neutron cross 
“Bsections of the common fissionable isotopes. 
further evidence is presented to show that U**® 
possesses the unusual negative-energy level 
required by the previous analysis. However, 
god fits are obtained for the cross sections of 
toth U5 and Pu?*® without such an unusual bound 
level, suggesting that the neutron resonance 
“fcross sections of the fissionable isotopes do not 
ahibita basic anomaly. The size of the level 
interference effects in each of the isotopes im- 
jlies that the fission process involves more 
tan one but no more than a few fission channels. 


YONTE CARLO CALCULATIONS OF NUCLEAR 
VAPORATION PROCESSES. IV. THE SPEC- 
TRAOF NEUTRONS AND CHARGED PARTICLES 
FROM NUCLEAR REACTIONS. Israel Dostrov- 
ay and Zeev Fraenkel, The Weizmann Institute 
Science, Rehovoth, Israel, and Lester Wins- 
‘erg, Lawrence Radiation Laboratory, Univer- 
tity of California, Berkeley, California, and 

The Weizmann Institute of Science, Rehovoth, 
B'szel (Received November 20, 1959). 


The calculation of spectra of neutrons and 
urged particles and of cross sections for their 
Mduction from nuclear reactions is compared 
th experimental values. A compound-nucleus 
wchanism followed by nuclear evaporation is 
“sumed for the reactions Zr, Ta, Bi(14.1-Mev 
iw); Ni(13.4-17.5-Mev n,p); Cu, Pd(23-Mev 
a); and Ni(162-Mev O*, a). The production of 


neutrons and charged particles from the interac- 
tion of 190-Mev protons with Ni, Ag, and Au is 
analyzed in terms of a nucleon cascade, followed 
by particle evaporation. The calculation of the 
nuclear evaporation is based on Weisskopf’s 
statistical theory. Fairly good agreement is ob- 
tained for the values of the cross sections for 
producing these particles with an appropriate set 
of radius and level-density parameters in each 
case. There are serious discrepancies, how- 
ever, in the comparison of the experimental and 
calculated spectra; many of the latter are defi- 
cient in low-energy neutrons and charged par- 
ticles. Possible improvements in the calcula- 
tion are discussed. 


MONTE CARLO CALCULATIONS OF NUCLEAR 
EVAPORATION PROCESSES. V. THE EMIS- 
SION OF PARTICLES HEAVIER THAN He‘. 

I. Dostrovsky, Z. Fraenkel, and P. Rabinowitz, 
The Weizmann Institute of Science, Rehovoth, 
Israel (Received November 20, 1959). 


Previous Monte-Carlo calculations of nuclear 
evaporation reactions have been extended to in- 
clude the emission of He®, Li®, Li’, Li®, and 
Be’ from Cu, Ag, Au, and Bi targets bombarded 
with high-energy protons (340-2000 Mev). Com- 
parison with available experimental results 
shows good agreement in most cases. A dis- 
crepancy has been observed between the calcu- 
lated and observed variation of Be’ formation 
cross section with the mass of the target nucleus, 
but even here the agreement is within a factor of 
two. It is shown that, for the usually chosen 
parameters of the calculation, a level density 
parameter of a=A/10 is necessary. 


TOTAL NEUTRON CROSS SECTIONS OF 
HELIUM, NEON, ARGON, KRYPTON, AND 
XENON. F. J. Vaughn, W. L. Imhof, R. G. 
Johnson, and M. Walt, Lockheed Missiles and 
Space Research Laboratory, Palo Alto, Cali- 
fornia (Received November 30, 1959). 


The total neutron cross sections of the noble 
gases helium, neon, argon, krypton, and xenon 
have been measured for neutron energies from 
120 kev to 6.2 Mev and from 12.1 Mev to 19.8 
Mev by a transmission experiment. The neu- 
trons were produced using the Li’(p,n)Be’, the 
T(p,n)He*, the D(d,n)He*, and the T(d, n)He* 
reactions in the appropriate energy intervals. 
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A Van de Graaff accelerator was the source of 
the protons or deuterons. In general, the re- 
sults obtained agree with previous work where 
such work exists. A previously unobserved 
S-wave scattering resonance was found in neon 
at about 500 kev, indicating the presence of an 
excited state in Ne* with J=1/2 and even parity. 
The results for argon, krypton, and xenon ex- 
hibit general agreement with the cross sections 
of neighboring elements, as would be expected 
from the previously observed smooth variation 
of the o(A, £) surface. 


HALF-LIFE OF Pb”°. W. R. Eckelmann,* W. S. 
Broecker, and J. L. Kulp, Lamont Geological 
Observatory, Columbia University, Palisades, 
New York (Received February 19, 1959). 


A new determination of the half-life of Pb*° 
has been made by the geological method. PbCl, 
extracted from uranium minerals in secular 
equilibrium was used to calibrate a thick-source 
scintillation counter for Po*’° alpha particles. 
Using this calibration, the absolute activity of 
Po”? in partial equilibrium with a known num- 
ber of Pb”° atoms prepared from the decay of 
a measured quantity of radon was determined. 
From these data a half-life of 21.4+0.5 years 
was obtained for Pb”"°. 


"Present address: Jersey Production Research 
Company, Tulsa, Oklahoma. 


DECAY OF Be®™ (2.43-Mev STATE). E. M. Hen- 
ley* and P. D. Kunz, Department of Physics, 
University of Washington, Seattle, Washington 
(Received July 13, 1959). 


The decay of the 2.43- Mev state of Be® is 
treated theoretically. Of the open two-body de- 
cay channels all but one involve a nuclear state, 
the energy of which is not well defined. The 
usual formalisms have been generalized to take 
this into account. The estimate of the decay 
rates is made by means of a variational internal 
wave function for the Be” state, based upon the 
alpha-particle model. It is found that the prin- 
cipal mode of decay is to He®+He*. Model- 
independent arguments are given to show that 
decay to the ground state of Be® should be inhib- 
ited. Furthermore, the momentum and angular 
distributions of alphas emitted in the decay 
through several two-particle decay modes are 
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computed. These latter calculations do not as. 
sume any specific nuclear model, but depend 
the weak assumption that the state is excited by 
a direct reaction. Comparison with recent 
measurements indicates that in addition to the 
He® - He* decay, approximately 7% of the decay 
occurs to the ground state of Be®, which is co. 
sistent with our calculations. 


* 
Presently on leave at the Institute for Theoretic 
Physics, Copenhagen, Denmark. 


CROSS SECTIONS FOR THE (n, 2m) REACTION 
IN N*, P*, Cu®, AND Pr“. J. M. Ferguson 
and W. E. Thompson, U. S. Naval Radiological 
Defense Laboratory, San Francisco, California 
(Received October 12, 1959). 


The (n,2n) cross sections have been measurei 
for N*“*, P**, Cu®, and Pr’, for neutron ener- 
gies from 12.5 to 18 Mev. The annihilation ra- 
diation emitted from the product nuclides was 
counted with two Nal(T1) crystals in coincidence, 
In the energy range measured, the cross sec- 
tions were found to vary, as follows: N", 3.03 
to 11.67 mb; P*, 0 to 74 mb; Cu, 186 to 836 
mb; Pr**!, 1231 to 1737 mb. The results are 
generally in agreement with those of others. 
The data are compared with curves plotted from 
Weisskopf’s theoretical expression for (n, 2n) 
cross sections. 


NUCLEAR ENERGY LEVELS OF Na™. Carl T. 
Hibdon, Argonne National Laboratory, Lemont, 
Illinois (Received November 5, 1959). 


The neutron cross-section data up to 350 kev 
show a number of relatively large peaks and 
many smaller ones among the 86 peaks observed, 
the widths ranging from 0.2 to 6 kev. Approxi- 
mately 50 small peaks were observed between 
60 and 200 kev. Above 200 kev, each of the pre- 
viously known peaks was resolved into two or 
more peaks and between these large peaks may 
narrower peaks were observed. The analyses 
show 9 s-wave levels and 46 p-wave levels, the 
remainder being d- and f-wave levels. A plot of 
the number of levels having energies <£,, 54 


function of the neutron energy E,, shows an essen: 


tially linear distribution of the levels. As 0b- 
tained from the reduced widths averaged over 
both values of J, the value of the strength func- 
tion for 7=0 is 0.06; averaged over all values of 
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Jfor 1=1 it is 0.65; and for higher values of / it 
js too large in comparison with the p-wave 


strength function. 





DIFFERENTIAL CROSS SECTIONS FOR NEU- 
TRON RESONANCE SCATTERING FROM Na”. 
R. 0. Lane and J. E. Monahan, Argonne National 
Laboratory, Lemont, Illinois (Received Novem- 
ber 12, 1959). 


The differential scattering cross section for 
neutrons on sodium has been measured with an 
energy spread of ~25 kev for neutron energies 
between 200 kev and 800 kev. The data are pre- 
sented in the form of Legendre polynomial 
coefficients. 


NEUTRON PRODUCTION BY HEAVY-ION BOM- 
BARDMENTS. Edward L. Hubbard, Robert M. 
Main,* and Robert V. Pyle, Lawrence Radiation 
Laboratory, University of California, Berkeley, 
California (Received September 11, 1959). 


Neutron yields from C™, N"*, and Ne”° bom- 
bardments of a number of target elements have 
been measured by an activation method. The 
maximum bombarding energies were 10.4 Mev 
per nucleon of the incident ion. Neutron yields 
have been calculated by assuming complete 
fusion of the two nuclei, with an interaction 
radius of y,=1.5x10"** cm, followed by de- 
excitation of the compound nuclei by neutron 
emission only. Calculated neutron yields are a 
factor of about two higher than experiment in 
the case of heavy target nuclei, with greater 
differences for light targets. Some possible re- 
finements of the theory that could bring the re- 
sults closer to agreement with experiment are 
mentioned. 


"Now at Tracerlab Inc., Richmond, California. 


PHOTONEUTRON CROSS SECTIONS OF Li, N, 
AND A. R. W. Fast, P. A. Flournoy, R. S. 
Tickle, and W. D. Whitehead, University of 
Virginia, Charlottesville, Virginia (Received 
November 13, 1959). 


Using a Halpern-type photoneutron detection 
system, the photoneutron yields from Li, N", 
and A*° have been measured as a function of the 








maximum bremsstrahlung energy from threshold 





to approximately 50 Mev. The method of Penfold 
and Leiss was used to extract from the yield 
curves the total neutron cross section: o7= 0(y,n) 
+ o(y, pn) + 2oly, 2n)+... . The results are com- 
pared with previous findings of other laboratories. 
No gross structure was detected in the lithium 
cross section in the giant resonance region. The 
data indicate that lithium has a high-energy tail 
on the cross section of considerable magnitude. 


TOTAL CROSS SECTIONS OF THE O'*(p, a)N*® 
AND 0'°(p,)F'® REACTIONS. J. M. Blair and 
J. J. Leigh, School of Physics, University of 
Minnesota, Minneapolis, Minnesota (Received 
November 9, 1959). 


The angular distribution of the a particles 
from the O"*(p, a)N** reaction and the total cross 
section for the 0'%(p,m)F*® reaction have been 
measured for proton energies between 2600 and 
3000 kev. Comparison of the total cross sec- 
tions for the two reactions eliminate the pre- 
viously observed differences in resonance ener- 
gies in these reactions. At most energies the 
a@ particles are preferentially emitted in the 
forward and backward directions and the angular 
distributions are rapid functions of proton en- 


ergy. 


ENERGY SPECTRA OF PROTONS FROM (4d, p) 
REACTIONS IN HEAVY ELEMENTS. B. L. 
Cohen, J. B. Mead, R. E. Price, K. S. Quisen- 
berry, and C. Martz, Radiation Laboratory, 
University of Pittsburgh, Pittsburgh, Pennsyl- 
vania (Received July 27, 1959). 


Surveys of proton energy distributions from 
(d, p) reactions were made on nuclei with Z >30 
using resolutions of 500 and 80 kev. The gross 
structure shows broad peaks due to the major 
nuclear shells, as expected from the fact that 
(d, p) stripping reactions excite single-particle 
states; peaks due to the subshell structure can 
be seen in some cases, especially in the heavier 
nuclei. The energies of the various peaks do 
not shift from element to element in the manner 
expected from simple theory; it is shown that 
this is not in conflict with neutron cross-section 
evidence, and possible explanations are pro- 
posed. The energy spacing between major shells 
derived from these measurements allows cal- 
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culation of the reduced mass for nucleons in 
nuclei; the result is very different from the 
predictions of Brueckner’s theory, but explana- 
tions for this discrepancy are advanced. The 
results and their interpretation given here are 
in direct conflict with the Wilkinson theory of 
gamma-ray giant resonances. The energy spec- 
tra are very similar at different angles, which 
indicates that the stripping process is predomi- 
nant at all angles. Deviations from Butler’s 
angular distribution theory at large angles must 
therefore be due to difficulties in that theory 
rather than due to the onset of a competing proc- 
ess. 


POLARIZATION OF 9-Mev PROTONS ELASTIC - 
ALLY SCATTERED FROM MAGNESIUM. A. B. 
Robbins, University of Birmingham, Birming- 
ham, England, and Rutgers University, New 
Brunswick, New Jersey, and G. W. Greenlees, 
University of Birmingham, Birmingham, England 
(Received November 30, 1959). 


The polarization as a function of angle has 
been measured for protons elastically scattered 
from a 1-Mev thick magnesium target with a 
mean energy of 9.1 Mev. The resulting polari- 
zation distribution is compared to a differential 
cross-section measurement with the same 
target. 


DEVIATIONS FROM THE AT =0 ISOTOPIC SPIN 
SELECTION RULE IN FERMI TRANSITIONS. 
Claude C. Bouchiat, Palmer Physical Laboratory, 
Princeton University, Princeton, New Jersey 
(Received September 23, 1959). 


Experimental deviations from the AT =0 iso- 
topic spin selection rule have been observed 
experimentally in J~J beta transitions. In the 
theory of a vector interaction with a conserved 
current these deviations have to be explained 
only in terms of isotopic spin impurities, while 
in the conventional theory exchange mesonic 
currents may also induce Fermi transitions with 
4ST #0. In this paper an attempt is made to esti- 
mate the contribution of the isotopic spin im- 
purities arising from the Coulomb interaction 
between the protons. The j -j coupling shell 
model is used to calculate the relevant Coulomb 
matrix elements. When all the nucleons outside 
the core are in the same orbit, the main con- 
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tribution comes from the Coulomb interaction high 
between the protons outside the core. A com. [inte 
parison between the empirical Fermi matrix mod 
element and the calculated one is performed in | ex¢! 
the case of Mn™, Sc**, Na™. The two quantities | m0 
agree fairly well for manganese-52. No such | ene! 
an agreement is found in the two other cases, | hig! 
the predicted M ; being too large for sodium-2 | fiss 
too small for scandium-44, at least by a factor /fiss 
ten. This discrepancy may reflect the inade- | er 
quacy of the j -j shell model to describe the . 
Coulomb effects or the presence of mesonic Stat 
effects. More experiments are needed to make 
a choice between these two possibilities. 

NU! 

NU 
DIRECT DETERMINATION OF INTERNAL Coy.| TU! 
VERSION COEFFICIENTS. Paresh Mukherjee, | Lat 
Saha Institute of Nuclear Physics, Calcutta, Net 
India (Received November 24, 1959). ma! 

Using the same source and instrument geome- | 1 
try, both the external as well as internal con- | ig 
version lines of the 1.17- and 1.33-Mev gamma | ° 
rays of Co® are scanned in a Siegbahn-Slatis [| ™¢ 
spectrometer. From the known internal con- { 
version coefficients of these gamma rays, the | 
instrument is calibrated for the direct deter- sin 
mination of the internal conversion coefficient } 
of any other gamma rays having energy near 1.3 | 
Mev. As an example, the internal conversion der 
coefficient of the 1.408-Mev gamma ray, in the {or 
decay of Eu’®*, is measured. The value obtained | © 
is 4.99x1074. - 

the 

ap 
FISSION OF Ra”*“* BY DEUTERONS AND HELIUM - 
IONS. R. C. Jensen* and A. W. Fairhall, Depart-f (, 
ments of Chemistry and Physics, University of dis 
Washington, Seattle, Washington (Received 
November 30, 1959). 

Fission induced in Ra*® by 14.5- and 22-Mev 
deuterons, and by 23.5-, 31-, and 43-Mev He M 
ions has been studied using radiochemical tech- } D 
niques. The mass distributions of fission B 
products for deuteron-induced fission is triple- Jy 
humped, corresponding to separate symmetric ve 
and asymmetric fission modes. The symmetric | ¢ 
mode dominates at the higher bombarding energy. | x 
The mass distributions observed for fission 
products from He ion-induced fission look more 
“normal”: asymmetric at the lowest bombard- . 





ing energy, becoming a single broad peak at the 
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highest bombarding energy. These results are 
interpreted in terms of a symmetric fission 
}mode which increases strongly with increasing 
excitation energy, and an asymmetric fission 
mode which occurs mainly at low excitation 
energies following neutron evaporation from 
highly excited compound nuclei. Asymmetric 
fission is interpreted to be disappearing as a 
)tission mode for nuclei of lower atomic num- 
ber than thorium. 








"Now at the Department of Chemistry, San Diego 
State College, San Diego, California. 


NUMERICAL EVALUATION OF THE PION- 
NUCLEON FORWARD SCATTERING AMPLI- 


_| TUDE. James W. Cronin, Palmer Physical 


Laboratory, Princeton University, Princeton, 
New Jersey (Received July 13, 1959; revised 
manuscript received January 7, 1960). 


The real part of the forward elastic scatter - 
ing amplitude for 1*-proton scattering has been 
evaluated from experimental cross sections by 
means of dispersion relations. Recent measure- 
ments indicate two peaks in the 7 -proton total 
cross section at 590 and 870 Mev incident pion 
kinetic energy. Tables of the real part of the 
forward scattering amplitude for 7*-proton 
scattering are presented as a function of inci- 
dent pion kinetic energy in the laboratory. The 
forward scattering amplitudes obtained from 
some recent 7+ scattering experiments are com- 
pared with the calculations. Measurement of 
the forward charge exchange cross section 
appears to be the most suitable way of investi- 
gating the predictions of the dispersion relation 
athigh energies. The possibility of detecting 
Coulomb interference at small angles is also 
discussed. 


NUCLEAR INTERACTIONS IN CARBON PRO- 
DUCED BY COSMIC RAYS WITH ENERGIES 
BETWEEN 10°° AND 10” ev. Luisa F. Hansen 
a W. B. Fretter, Department of Physics, Uni- 
versity of California, Berkeley, California (Re- 
teived November 16, 1959; revised manuscript 
received February 3, 1960). 


An experiment is described in which high- 
‘ergy nuclear interactions in the range of 
*ergies 10*°- 10'? ev were analyzed by means 












of a cloud chamber in a magnetic field. Meas- 
urements of ionization and momentum made 
possible the identification of electrons and 7 
mesons to about 20 Bev/c. Protons, K mesons, 
and hyperons could not be identified unambigu- 
ously among themselves, except in very limited 
regions of momentum. The primary particles 
were cosmic-ray nucleons and a possible frac- 
tion of pions, the target nuclei were carbon, and 
the velocities of the primaries were determined 
from balance of momentum in the center-of-mass 
system. A total of 41 events were analyzed, and 
the results compared to previous experimental 
work and the predictions of the theories of Heisen- 
berg and Landau. The measurements made in- 
cluded the transverse momenta of the secondar - 
ies and their average energy in the center -of- 
mass system, the energy and angular distribu- 
tion of the pions and heavy particles (protons, 

K mesons, hyperons) in the center-of-mass sys- 
tem, the inelasticity of the collision, the mul- 
tiplicity of the showers, the percentage of 
strange particles, and the positive excess of the 
secondaries. 


HIGH -ENERGY K~-MESON INTERACTIONS AND 
DECAYS. Stanley C. Freden, Francis C. Gilbert, 
and R. Stephen White, Lawrence Radiation Labor - 
atory, University of California, Livermore, 
California (Received October 20, 1959). 


The interactions of 20- to 300-Mev K~ mesons 
on free protons and on emulsion nuclei have been 
studied and their decays analyzed. The K~-meson 
mean lifetime is found to be (1.38 + 0.24) x10-8 
sec. Examples of the decay modes Ky2, Ky2; 
Ky3, 7> and 7’ are identified. The branching 
ratios are found to be in agreement with those 
for K* mesons. The (K~p,K~p) elastic scatter - 
ing cross section is found to be 35+ 16 mb and 
the (K~p, £*2*) cross section 27+13 mb in the 
energy region of 150 to 300 Mev. The mean 
free path for K~-meson captures and inelastic 
scatters on emulsion nuclei, except hydrogen, 
can be represented by A(cm) = (17.2 +3.4) 
+(0.081+0.027)T,, for K~-meson energies from 
20 to 300 Mev. This increase in the mean free 
path with K~-meson energy is explained in terms 
of the decreasing nucleon cross section. The 
fraction of the interactions of K~ mesons on 
nuclei which are inelastic scatters increases 
from 2% at low K~-meson energies to about 
15% at 150 Mev. This increase in inelastic scat- 
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tering with energy is additional evidence that the 


nucleus is becoming partially transparent to 


K~ mesons at about 150 Mev. Data are presented 
for the fraction of events with observed 7 mesons, 
= hyperons, and 7 meson-z hyperon pairs. The 
data are discussed in terms of the model which 


was previously presented to explain K"-meson 
captures at rest on nucleons bound in nuclei. 


a* -p INTERACTIONS PRODUCING =* - K* IN 
A PROPANE BUBBLE CHAMBER. W. H. 
Hannum, * H. Courant, E. C. Fowler, H. L. 
Kraybill, J. Sandweiss, and J. Sanford, Yale 


University, New Haven, Connecticut (Received 


November 13, 1959). 


Seven cases of production of £* - K* by 1.0- 
Bev n° - p interactions have been identified in 
the Yale propane bubble chamber. In six of 
these cases, the =* direction was forward in 
the center-of-momentum system. The cross 


section derived from these is 0.08 mb. The de- 
cay of the =* is equally divided between the pion 


and proton decay modes. 


*Now at Westinghouse Atomic Power Laboratory, 
Pittsburgh, Pennsylvania. 


FURTHER SEARCH FOR THE DECAY i*— er +y. 
S. Frankel, V. Hagopian, J. Halpern, and A. L. 


Whetstone, University of Pennsylvania, Phila- 
delphia, Pennsylvania (Received November 4, 
1959). 


A new experiment for determining the upper 
limit for the branching ratio R of the process 
ut~—e*++y relative to the normal decay mode 
u*t~—e*+v+PD yields a value of R of less than 
1.2x10-*® with a 90% confidence level. 


EVIDENCE FOR A FORBUSH TYPE OF DE- 
CREASE IN THE INTENSITY OF HEAVY NU- 
CLEI OF THE PRIMARY COSMIC RADIATION. 


S. Biswas,* P. J. Lavakare,f K. A. Neelakantan, 


and P. G. Shukla, Tata Institute of Fundamental 
Research, Bombay, India (Received November 
13, 1959). 


In an emulsion stack flown on March 13, 1956, 


from Iowa, the flux of heavy nuclei with Z 26 in 
the primary cosmic radiation was measured 
as 15.8 and 4.7 particles/m? sec sr for par- 
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ticles of kinetic energy 2 0.23 and > 1.55 Bev/ 
nucleon, respectively. The measured flux of 
energy > 1.55 Bev/nucleon was 45+ 11% lower 
than the normal flux. It is shown that almost th 
entire part of the reduction must be attributed 


to a large Forbush decrease of the cosmic rag. 


tion that occurred at the same time. The ex- 
ponent of the integral energy spectrum of heavy 
nuclei (Z > 6) was measured as 1.74 0.25 in the 
energy interval 0.23 to 9 Bev/nucleon. As this 
value is not significantly different from its nor. 
mally measured value, it appears that the large 
reduction in the primary flux was not accompa- 
nied by any significant change in the energy 
spectrum. 


“Now at the University of Minnesota, Minneapolis, 
Minnesota. 

tNow at the University of Rochester, Rochester, 
New York. 


ANGULAR DISTRIBUTION OF FRAGMENTS 
FROM THE FISSION OF BISMUTH BY 450- Mev 
PROTONS. M. V. Ramaniah* and Nathan Sugar- 
man, Enrico Fermi Institute for Nuclear Studies, 
University of Chicago, Chicago, Illinois (Re- 
ceived November 4, 1959). 


Test experiments on the “bead-cone” method 
used for the measurement of the angular distri- 
bution of fission fragments from the 450-Mev 
proton fission of bismuth have been performed 
These experiments make doubtful the conclusion 
of the previous measurement. 


* 

Present address: Radiochemistry Laboratory, 
Atomic Energy Establishment, Trombay, Bombay, 
India. 


SMALL-ANGLE PROTON SCATTERING AT 3 
Bev. W. M. Preston, Richard Wilson, and J. C. 
Street, Cyclotron Laboratory, Harvard Univer- 
sity, Cambridge, Massachusetts (Received 
November 2, 1959). 


The differential cross section for elastic scat- 
tering of 3-Bev protons has been measured from 
targets of hydrogen, carbon, copper, and lead 
over the angular range 0.5 to 4 degrees in the 
laboratory coordinate system. Within our limits 
of error, no evidence of Coulomb-nuclear inter 
ference was found. In the case of hydrogen, the 
real component of the nuclear scattering ampli- 
tude (if any) is not over one-tenth the imaginary; 
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less stringent limits can be set for the other 
elements because of experimental uncertainties. 
For carbon, the differential nuclear scattering 
cross section at zero degrees agrees with that 
predicted by the optical theorm for a purely 
imaginary scattering amplitude; the lack of such 
agreement found for hydrogen may be evidence 
for a spin-dependent term. A derived value of 
the rms radius for p-p scattering exceeds that 
found in electron-proton scattering by a factor 


12. 


INTERACTION OF K*+ MESONS WITH PROTONS. 
T. F. Kycia,* L. T. Kerth, and R. G. Baender, 
Lawrence Radiation Laboratory, University of 
California, Berkeley, California (Received 
September 11, 1959). 


The total K* - p cross section was measured at 
the three K*-meson energies 175425, 225425, 
and 275+25 Mev, and the differential scattering 
cross section was measured at 225 Mev. The 
k*-p nuclear force was shown to be repulsive, 
from the observed constructive interference with 
Coulomb scattering. The differential cross sec- 
tion was otherwise isotropic and could arise 
from either pure S-wave or pure P-wave scat- 
tering. 

Subtracted dispersion relations were applied 
to these data and the rest of the available K- 
proton scattering data. The statistical errors 
inthe data were found to be too large to deter - 
mine the K-hyperon relative parity. On the 
assumption that the KA and KZ relative parities 
are the same, then for scalar coupling, g*/47 
would be less than 0.6; for pseudoscalar coup- 
ling, it would be less than 10. 


"Present address: Brookhaven National Laboratory, 
Upton, New York. 


NEUTRON FORM FACTORS FROM HIGH- 
ENERGY INELASTIC ELECTRON-DEUTERON 
SCATTERING. S. Sobottka,* Department of 
Physics and High-Energy Physics Laboratory, 
Sanford University, Stanford, California (Re- 
telived November 30, 1959). 


The inelastic electron-deuteron scattering 
‘toss section has been measured for incident 
‘lectron energies between 300 Mev and 650 Mev 
and for final electron energies primarily at the 
taxima of the inelastic continuua. The data 





were interpreted in terms of neutron form fac- 
tors by employing the impulse approximation 
calculations of Goldberg. The results indicate 
that Fo,” is nearly equal to the proton form 
factor Fp’ for 2.65<q?<15.1 (fermi)~* but may 
be 20% or 30% higher than Fp° for the lowest 
of these g values. Uncertainties, primarily in 
the theory, make it impossible to determine 
whether the difference is real. The results also 
indicate that -2.5<F4,/F9, <0.5 for 5.1<q?<12.8 
e~, 

"Now at Boeing Scientific Research Laboratories, 
Boeing Aircraft Company, Seattle, Washington. 


COSMIC -RAY HEAVY NUCLEI AT THE GEO- 
MAGNETIC EQUATOR. D. D. Kerlee and O. K. 
Krienke, Jr., Seattle Pacific College Institute 

for Research, Seattle, Washington, J. J. Lord, 
Department of Physics, University of Washing- 
ton, Seattle, Washington, and M. E. Nelson, 
Department of Physics, College of Puget Sound, 
Tacoma, Washington (Received October 16, 1959). 


The flux of primary cosmic-ray heavy nuclei 
has been measured with a sandwich of C-2 and 
G-5 emulsions near the island of Guam, mag- 
netic latitude 4°N. The emulsions were exposed 
at about 101 000 ft for 7.3 hours and the calculated 
fluxes at the top of the atmosphere in particles/ 
m*-sec-sr were found to be: 


Li, Be,B (L nuclei) 0.26+0.05, 
C,N,O,F (M nuclei) 0.66+0.10, 
Z210 (HW nuclei) 0.37+0.09. 


This flight took place on February 12, 1957, 
about a year prior to the sun spot maximum; 
however, there is no evidence within statistical 
error for any changes in the composition of the 
primary heavy nuclei. 


HIGH ORBITAL S-STATE CAPTURE OF 17 
MESONS BY PROTONS. T. B. Day, G. A. Snow, 
and J. Sucher, Department of Physics, Univer- 
sity of Maryland, College Park, Maryland (Re- 
ceived November 30, 1959). 


The consequences of the very short capture 
time for 7 mesons in liquid hydrogen, recently 
measured by Fields, Yodh, Derrick, and Fetko- 
vich, are discussed. It is pointed out that colli- 
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sional de-excitation mechanisms, even including 
the Stark effect enhancement of capture, are in- 
adequate to explain the experiment. Alternative 
possibilities are discussed. 


MUON CAPTURE IN He’. Akihiko Fujii, Depart- 
ment of Physics, Purdue University, Lafayette, 
Indiana (Received November 20, 1959). 


The hard-core wave function for a three- 
nucleon system is used to calculate the capture 


rate of the reaction » ~ +He* ~ H°(ground state) + v. 


It is found to be 1.66x10° sec™. 


INTERFERENCE PHENOMENA IN NUCLEAR 
SCATTERING OF NEUTRAL K MESONS. 
Nripendra N. Biswas, Lawrence Radiation La- 
boratory, University of California, Berkeley, 
California (Received November 16, 1959). 


The scattering of neutral K mesons has been 
treated phenomenologically. The scattered 
beam, in general, contains both K, and K, com- 
ponents having different amplitudes. These 
amplitudes interfere with each other in the 
generation of K, and K, components in the scat- 
tered beam. The relative sign of the two ampli- 
tudes may then be determined from the analysis 
of K,, K, decays. The leptonic decay rates of 
the scattered beam show a dependence on AM, 
the mass difference between K, and K,, in such 
a way that the sign of 4M can, in principle, be 
determined experimentally. 


LIFETIME MATRIX IN COLLISION THEORY. 
Felix T. Smith, Stanford Research Institute, 
Menlo Park, California (Received October 16, 
1959). 


The duration of a collision is usually a rather 
ill-defined concept, depending on a more or less 
arbitrary choice of a collision distance. If the 
collision lifetime is defined as the limit, as 
R-o, of the difference between the time the 
particles spend within a distance R of each other 
and the time they would have spent there in the 
absence of the interaction, a well-defined quan- 
tity emerges which is finite as long as the inter- 
action vanishes rapidly enough at large R. 

In quantum mechanics, using steady-state 
wave functions, the average time of residence 
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in a region is the integrated density divided by 
the total flux in (or out), and the lifetime is 
defined as the difference between these resideng: 
times with and without interaction. Transform. 
tion properties require construction of the life- 
time matrix, Q. If the wave functions 7; are 
normalized to unit total flux in and out through 
a sphere at R~o, the matrix elements are 


li " *dt -R(v.76 
ae oe 0,0, dr -Rlo, 6, 


” * 

2S ine Six | 
av 

where the average value is taken to eliminate 
oscillating terms at large R. Here, Sjp is an 
element of the unitary scattering matrix, S, 
and v; is the velocity in the ith channel. Q is 
Hermitian; a diagonal element Q;; is the aver- 
age lifetime of a collision beginning in the ith 
channel. As a function of the energy, Q is 
related to S: Q= -ihSdStdE; Q and S contain the 
same information, from complementary points 
of view. When Q is diagonalized, its proper 
values, q;;, are the lifetimes of metastable 
states if they are large compared to #/E; fora 
sharp resonance, the measured lifetime is the 
average of g;;(E) over a distribution in energy. 
The corresponding eigenfunctions, W;, are the 
proper functions to describe these metastable 
states. The causality principle appears directly 
from an inequality involving the integral expres- 
sion for Qj; or 9;;, and it is shown how some 
of its consequences for inelastic collisions can 
be deduced. 


RADIATIVE PION DECAY INTO ELECTRONS. 
Sidney A. Bludman and James A. Young, Law- 
rence Radiation Laboratory, University of Cali- 
fornia, Berkeley, California (Received June 25, 
1959; revised manuscript received January 8, 
1960). 


The possibility of distinguishing the pion struc 
ture- dependent radiation from the conventional 
inner bremsstrahlung radiation in the radiative 
decay of pions into electrons is discussed. Cal- 
culation of the photon energy spectrum and 
angular correlation shows that evidence for pio 
structure would be obtained if any photons of el 
ergy less that 70 Mev were detected in 180° 














coincidence with s-decay electrons. The proba- 
bility of such events per unit solid angle is 20.2 
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x10"? relative to ordinary 1— 1+ v decay, if the 
assumption of a conserved vector current is 

made to relate the rate of radiative decay through 
the weak V interaction to the rate of 1°~ 2, decay. 


THEORY OF ALLOWED AND FORBIDDEN 
TRANSITIONS IN MUON CAPTURE REACTIONS. 
Masato Morita, Columbia University, New York, 
New York, and Akihiko Fujii, Brookhaven 
National Laboratory, Upton, New York (Re- 
ceived November 9, 1959). 


A general formula for the transition rate of 
the muon capture reaction, u~+(A,Z)—v+(A,Z-1), 
where the final nuclear state has definite spin and 
parity, is given in terms of the total and orbital 
angular momenta of the emitted neutrino and of 
the spins of the initial and final nuclear states. 
The induced pseudoscalar interaction and the in- 
teraction due to the assumption of conserved 
vector current are taken into account, together 
with the vector and axial vector interactions. 
The forbiddenness of the muon capture reaction 
is defined in a manner analogous to the theory 
of the beta decay. The spin and parity changes 
can assume the values (0+, 1+), (0-,1-,2-), 
In(-/", n +1(-)"] for the allowed, first forbidden, 
and nth (mn 22) forbidden transitions, respectively. 
(+and - mean the parity change “no” and “yes.”) 
For these transitions, the number of reduced 
nuclear matrix elements involved is nine, six- 
teen, and fourteen, respectively. The transition 
rate of the muon capture reaction is reduced by a 
factor of 10°, approximately, for a two-unit in- 
crease of the forbiddenness, if the atomic num- 
ber and the energy of the neutrino are constant. 
The contribution from the higher order transi- 
tion to the lower one is less than 0.1% in the 
medium and light nuclei. Explicit formulas for 
the transition rate are given for the allowed, 
first forbidden, and nth forbidden transitions. 
They are related to the corresponding formulas 
of beta decay. Our formalism was applied to the 
calculation of the partial muon capture rate by 





C™ ending in the ground state of B’?. The numer- 
ical analysis indicates that measurement of this 
capture rate can determine whether the con- 
served vector current interaction term exists 
in nature only if the coupling constant of the in- 
duced pseudoscalar interaction and the nuclear 
wave functions are well known. The transition 
rates are given for the j-j coupling shell model 
and harmonic oscillator wave functions. They 
are 9-13% smaller than those given by Fujii and 
Primakofi. 


ELECTROMAGNETIC STRUCTURE OF THE 
NUCLEON. Adam M. Bincer, Department of 
Physics, University of California, Berkeley, 
California (Received November 16, 1959). 


Dispersion relations are proved for the elec- 
tromagnetic and mesonic nucleon vertex func- 
tions considered as a function of the nucleon 
mass. The results are used to express the 
isotopic scalar and the isotopic vector electro- 
magnetic form factors of the nucleon in terms 
of pion electroproduction (or photoproduction) 
and pion-nucleon scattering amplitudes in the 
J=2,7T =3 state. 
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T, P, C SYMMETRIES IN THE 7° DECAY. 
Jeremy Bernstein and Louis Michel,* Ecole 
Polytechnique, Paris, France (Received Novem- 
ber 19, 1959). 


An analysis is given of the decay of the 7° in 
which allowance is made for possible break- 
downs in T, P, and C symmetries. It is shown 
that experiments, until now, have demonstrated 
only that the two-photon state is an eigenstate 
of 7P, but not of T and P separately. A dis- 
cussion of experiments which may verify T and 
P symmetry for the two-photon state is given. 


*Present address: Service de Physique Théorique, 
Faculté des Sciences, B. P., Orsay (Seine-et-Oise), 
France. 





